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»' ' A theory is formulated for time dependent fluctuations of the spectrum of a single molecule in 

a dynamic environment. In particular, we investigate the photon counting statistics of a single 
^ ' molecule undergoing a spectral diffusion process. Based on the stochastic optical Bloch equation, 

Q I fluctuations are characterized by Mandel's Q parameter yielding the variance of number of emitted 

photons and the second order intensity correlation function, (ji'^-'(f). Using a semi-classical approach 
and linear response theory, we show that the Q parameter can be described by a three-time dipole 
0^ ' correlation function. This approach generalizes the Wiener-Khintchine formula that gives the aver- 

age number of fluorescent photons in terms of a one-time dipole correlation function. We classify 
the time ordering properties of the three-time dipole correlation function, and show that it can be 
represented by three different pulse shape functions similar to those used in the context of nonlinear 
^ , spectroscopy. An exact solution is found for a single molecule whose absorption frequency under- 

. goes a two state random telegraph process (i.e., the Kubo-Anderson sudden jump process.) Simple 

expressions are obtained from the exact solution in the slow and fast modulation regimes based on 
appropriate approximations for each case. In the slow modulation regime Q can be large even in 
the long time limit, while in the fast modulation regime it becomes small. 



I. INTRODUCTION 



> 

^ • In recent years, a new approach to condensed phase spectroscopyhas emerged, that focuses on the spectral properties 
of a single molecule (SM) embedded in a condensed phaseP,|2,ll,0|- Thanks to experimental advances made in optics 

■ and microscopy's! , it is now possible to perform single molecule spectroscopy (SMS) in many different systems. 
O I Motivations for SMS arise from a fundamental point of view (e.g., the investigation of the field-matter interaction 

at the level of a SM, the verification of statistical assumptions made in ensemble spectroscopy, etc) and from the 
possibility of applications (e.g., the use of SMS as a probe for large biomolecules for which a SM is attached as a 
y—i fluorescent marker). 

^ In general, the spectral properties of each individual molecule vary from molecule to molecule due to differences in 

00 the local environments with which each SM is interacting[l[l|l|l[l3,[l3,[ll|ll[ll[ll[ll. With its unique ability to 

■ detect dynamical phenomena occurring at the level of an individual molecule surrounded by its local environment, SMS 
I ] has uncovered the statistical distributions of microscopic quantities of the environment that are hidden in traditional 

. ensemble averaged spectroscopy. In particular, a single molecule spectrum measured for a finite time necessarily "sees" 
' the temporal fluctuations of the host environment that occur on timescales comparable to the measurement timescale, 
. and therefore lead, in many cases, to a stochastically fluctuating siiigle molecule spectrum. Time dependent fluctuation 
' phenomena in SMS occur in many ways, such as spectral diffusion p, 1^ ISi, .9,, ,1 2,, _ 1 7j and fluorescence intermittencvflM 
J3 ■ ESi 0j nil ■ The physical mechanisms causing these fluctuation phenomena vary depending on the dynamical 

processes a SM is underg oing , including: triplet state dvnamics llSl Il9|. l2Ci . energy transfer processes l24l |25|. [26l . 
exciton transfer proces ses Il7l ! chemical reactions l27l I2M f29l Isol l3lll32j. conformational changes S^, l34l ISSl l36l |. 
rotational dvnamics|37ll38l|. and diffusion processes |39|. Thus SMS provides a unique microscopic tool to investigate 
Q ' the dynamical processes that a SM and its environment undergo during the measurement time. 
^ I One important process responsible for time-dependent fluctuations in SMS is spectral diffusion, i.e., perturbations 
' or excitations in the environment of the SM produce random changes in the transition frequency of the SM 1^,0,0, 
• ^ . M, Q> ^> ^> nil 1^ J leading to a time-dependent spectrum. Spectral diffusion processes have been observed 
^ ■ in various systems including dye molecules in a molecular crystal00| and in a low temperature glassH,0], quantum 



H ' dotsl45|, light harvesting svstems |46ll47| . and dendrimers'48|. Since the spectral diffusion process directly reflects both 
(i) the interaction between the SM and its environment and (ii) the local dynamics of the latter, careful analysis of the 
time-dependent fluctuations of SMS illuminates the interplay between various dynamical processes in the condensed 
phase. In this work we formulate a stochastic theory of SMS undergoing a spectral diffusion process. In particular, 
we address the issue of the counting statistics of emitted photons produced by a SM undergoing a spectral diffusion 
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process. Our studies show how the fluctuations in SMS can be used to probe the dynamics of SM and its interaction 
with the excitations of the environment. 

Previously, the photon counting statistics of an ensemble of molecules, studied by various methods, for example, 
the fluorescence correlation spectroscopy^^ls^, 1^115, '53] has proved useful for investigating dynamical processes of 
various systems. The photon statistics of a SM is clearly different from that of the ensemble of molecules due both to 
the absence of inhomogcneous broadening and to the correlation between fluorescence photons that exists only on the 
SM level. In some SMS experiments, the measurement time is limited due to photobleaching, where the emission of a 
SM is quenched suddenly because of various reasons, for example, reaction with oxygen. Thus it is not an easy task 
in general to collect a sufficient number of photon counts to have good statistics. However, many SM-host systems 
have been found to remain stable for long enough time to measure photon statistics [s^ 0, IsMIs^ Is^ |5^ [s^. In 
view of these recent experimental activities, a theoretical investigation of the counting statistics of photons produced 
by single molecules, in particular when there is a spectral diffusion process is timely and important. 

We will analyze the SM spectra and their fluctuations semiclassically using the stochastic Bloch equation in the 
limit of a weak laser field. The Kubo-Anderson sudden jump approach 60), 61,, 62!| is used to describe the spectral 
diffusion process. For several decades this model has been a useful tool for understanding line shape phenomena, 
namely, of the average number of counts (n) per measurement time T, and has found many applications mostly in 
ensemble measurements, e.g., NMR,|63| and nonlinear spectroscopy. ,64j More recently, it was app lied to model SMS 
in low temperature glass systems in order to describe the static properties of line shapes [T^ Il5l ITgL l65| and also to 
model the time dependent fluctuations of SMS. 66, 67, 68] 

Mandel's Q parameter quantitatively describes the deviation of the photon statistics from the Poissonian case [6^[70l] . 



Q= "''7;'" -1, (1) 

(n) 

where n is the random number of photon counts, and the average is taken over stochastic processes involved. In the 
case of Poisson counting statistics (5 = while our semiclassical results show super-Poissonian behavior (Q > 0) for 
a SM undergoing a spectral diffusion process. For short enough times, fluorescent photons emitted by a SM show 
anti-bunching phenomena (— 1 < Q < 0), a sub-Poissonian nonclassical effect ,,55» .56.. JI. ..72^ .73„ .74] . Our semiclassical 
approach is valid when the number of photon counts is large. Further discussion of the validity of our approach is 
given in Section IVlIII 

One of the other useful quantities to characterize dynamical processes in SMS is the fluorescence intensity correlation 
fimction, also called the second-order correlation fimction, g^^\t), defined bv fTsLlT^ 

.'"W = ^ffi^- (2) 

This correlation function has been used to analyze dynamical processes involved in many SMS experiments!^ |2^ 
E^lzEIzS' Here I{t) is the random fluorescence intensity observed at time t. It is well known that for a stationary 
process there is a simple relation between g'-^^t) and Q(T) 55. 73..80j 

Q{T) -^1^ dii dt2.g(2) it2) - {I)T, (3) 

where T is the measurement time. 

The essential quantity in the present formulation is a three-time correlation function, C3 (ri, T2,_T3), which is similar 
to the nonlinear response function investigated in the context of four wave mixing processes |64|. The three-time 
correlation function contains all the microscopic information relevant for the calculation of the lineshape fluctuations 
described by Q. It has appeared as well in a recent paper of PlakhotnikjH3 in the context of intensity-time-frequency- 
correlation technique. In the present work, important time ordering properties of this function are fully investigated, 
and an analytical expression for Q is found. The relation between C3 (ti, T2,t^) and lineshape fluctuations described 
by Q generalizes the Wiener-Khintchinc theorem, that gives the relation between the one-time correlation function 
and the averaged lineshape. 

The timescale of the bath fluctuations is an important issue in SMS. Bath fluctuations are typically characterized 
as being in either fast or slow modulation regimes (to be defined later) ^65] . If the bath is very slow a simple adiabatic 
approximation is made based on the steady state solution of time-independent Bloch equation. Several studies have 
considered this simple limit in the context of SMSjMl^ fifiUZll ]. From a theoretical and also experimental point of view 
it is interesting to go beyond the slow modulation case. In the fast modulation case it is shown that a factorization 
approximation for the three-time correlation function yields a simple limiting solution. In this limit the lineshape 
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exhibits the well known behavior of motional narrowing (as timescale of the bath becomes short, the line is narrowed). 
By considering a simple spectral diffusion process, we show that Q exhibits a more complicated behavior than the 
lineshape does. When the timescale of the bath dynamics goes to zero, we find Poissonian photon statistics. Our 
exact results can be evaluated for an arbitrary timescale of the bath and are shown to interpolate between the fast 
and slow modulation regimes. 

This paper is organized as follows. In Sec. Ill Al the stochastic Bloch equation is presented and a brief discussion 
of its physical interpretation is given, and in Sec. Ill Bl the prescription for the relation between the solution of the 
optical Bloch equation and the discrete photon counts is described. We briefly review several results on counting 
statistics, which will later clarify the meaning of some of our results. Section Hill presents simple simulation results of 
SM spectra in the presence of the spectral diffusion to demonstrate a generic physical situation to which the present 
theory is applicable. In Sec. lIVI an important relationship between Q and the three-time correlation function is found, 
and the general properties of the latter are investigated. An exact solution for a simple spectral diffusion process is 
found in Sec.0 In Sec. I VII we analyze the exact solution in various limiting cases so that the physical meaning of our 
results becomes clear. Connection of the present theory to experiments is made in Sec. IVIll In Sec. IVIIll we further 
discuss the validity of the present model in connection with other approaches. We conclude in Sec. IIXI Many of the 
mathematical derivations are relegated to the Appendices. 



Our theory presented in this section consists of two parts; first, we model the time evolution of a SM in a dynamic 
environment by the stochastic optical Bloch equation, and second, we introduce the photon counting statistics of a 
SM by considering the shot noise process due to the discreteness of photons. 



We assume a simple nondegenerate two level SM in an external classical laser field. The electronic excited state |e) 
is located at energy cjq above the ground state \g). We consider the time-dependent SM Hamiltonian 



where cr^ is the Pauli matrix. The second term reflects the effect of the fluctuation of the environment on the absorption 
frequency of the SM coupled to J perturbers. The stochastic frequency shifts ^(^j{t) (i. e. the spectral diffusion) are 
random functions whose properties will be specified later. The last term in Eq. describes the interaction between 
the SM and the laser field (frequency io^), where d = d^^gax is the dipole operator with the real matrix element 
deg ~ {e\d\g). We assume that the molecule does not have any permanent dipole moment either in the ground or in 
the excited state, (gidig) = (e|d|e) = 0. 

In the limit of a weak external field the model Hamiltonian describes the well known Kubo-Anderson random 
frequency modulation process whose properties are specified by statistics of Aujj{t) 60, 61, 62]. When the fluctuating 
part of the optical frequency Aloj is a two state random telegraph process, the Hamiltonian describes a SM (or spin 
of type A) coupled to J bath molecules(or spins of type B), these being two level systems. Under certain conditions 
this Hamiltonian describes a SM interacting with many two level systems in low temperature glasses that has been 
used to analyze SM lineshapcs 14, 15, 16, 65, 67, ..681] . 

The molecule is described by 2 x 2 density matrix p whose elements are pgg , pee , Pge , and peg ■ Let us define 



II. THEORY 



A. Stochastic Optical Bloch Equation 




(4) 



u = 



V — 



-(pgee-^--* + 
— (n p"''^^'* - 



Pege' 



) 



) 



(5) 



(6) 



w = 



„ (Pee Pgg): 



(7) 



and note that from the normalization condition pee + Pgg = 1, we have pee = w + 1/2. By using Eq. the stochastic 
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Bloch equations in the rotating wave approximation are given bv |8lL l8 

u ^ 5L{t)v-^-^, (8) 
Vv 

V = -^5L{t)u - — -flw, (9) 

r 

w = Qv-Tw- -. (10) 

1/r is the radiative lifetime of the molecule added phenomenologically to describe spontaneous emission, O — —dgg ■ 
Eq/H is the Rabi frequency, and the detuning frequency is defined by 

SL{t) = LUL-L^o- ^io{t), (11) 

.7 

Au;{t) - ^A^,(i). (12) 

i=i 

Besides the natural relaxation process described by F, other Ti and T2 processes can easily be included in the present 
theory, w represents half the difference between the populations of the state |e) and \g), while u and v give the mean 
value of the dipole moment d, 

Tr (pd) — 2dge [u cos (wl^) — v sin {t^Lt)] ■ (13) 

In recent studies [s^ Is^ it has been demonstrated that the deterministic two level optical Bloch equation approach 
captures the essential features of SMS in condensed phases, which further justifies our assumptions. 

The physical interpretation of the optical Bloch equation in the absence of time-dependent fluctuations is well 
known0, [t^- Now that the stochastic fluctuations are included in our theory we briefly discuss the additional 
assumptions needed for standard interpretation to hold. The time-dependent power absorbed by the SM due to work 
of the driving field is, 

dW d 

— =cosiujLt)Eo--TT{pd). (14) 
at at 

As usual, additional averaging (denoted with overbar) of Eq. (|14|l over the optical period of the laser is made. This 
averaging process is clearly justified for an ensemble of molecules each being out of phase. For a SM, such an 
additional averaging is meaningful when the laser timescale, 1/lul, is much shorter than any other timescale in the 
problem (besides I/ujq, of course). By using Eq. H13(l this means. 



VCOs'^{uJLt) ~ VCOS^{u!Lt)j (15) 

under the conditions \u{t)\ <C ijJL|w(t)| etc, and hence we have 

^ ~ nnujLvit). (16) 

The absorption photon current (unit l/[time]) isl75| 

1 dW 

m = —-^ = nv{t). (17) 
niOL at 

Neglecting photon shot noise (soon to be considered), J^^ I{t)dt has the meaning of the number of absorbed photons 
in the time interval (0,T) (i.e., since Jq {^^)dt is the total work and each photon carries energy huL). By using 
Eqs. (jHl)-(|inil, we have 

Pee = - Fpee- (18) 

In the steady state, pee — 0, we have flv — Tpee, and since Slu has a meaning of absorbed photon current, Tpee has 
the meaning of photon emission current. For the stochastic Bloch equation, a steady photon flux is never reached; 
however, integrating Eq. H18fl over the counting time interval T, 



[PeeiT) - peem +T ( Pee{t')dt' = 17 / v{t')dt' , 

Jo Jo 



(19) 



5 



and using \pf,e{T) — pee(0)| < 1 we find for large T that the absorption and emission photon counts are approximately 
equal, 

r ( p,,{t)At ~ / I[t)At, (20) 
Jo Jo 



emitted photons absorbed photons 

provided that I{t)dt ^ 1. Eq. (|20|l is a necessary condition for the present theory to hold, and it means that 
the large number of absorbed photons is approximately equal to the large number of emitted photons (i. e. we have 
neglected any non-radiative decay channels). When there are non-radiative decay channels involved, one may modify 
Eq. (|20|) approximately by taking into account the fluorescence quantum yield, 0, the ratio of the number of emitted 
photons to the number of absorbed photons, 

Pee(t)dt~0 / I{t)dt. (21) 



B. Classical Shot Noise 



Time dependent fluctuations are produced not only by the fluctuating environment in SMS. In addition, an impor- 
tant source of fluctuations is the discreteness of the photon, i.e., shot noise. Assuming a classical photon emission 
process, the probability of having a single photon emission event in time interval (t, t + dt) isj70| 

Prob(t, t + dt) ^ rpee{t)dt. (22) 

While this equation is certainly valid for ensemble of molecules all subjected to a hypothetical identical time-dependent 
environment, the validity of this equation for a SM is far from being obvious. In fact, as we discuss below, only under 
certain conditions we can expect this equation to be valid. By using Eq. H22|l the probability of recording n photons 
in time interval (0, T) is given by the classical counting formula 70] 

W" 

p{n,T) = —exp{-W), (23) 
n\ 

with 

W = rj [ I{t)dt, (24) 
Jo 

where 77 is a suitable constant depending on the detection efficiency. For simplicity we set rj — I here, but will 
re- introduce it later in our final expressions. Here W = rjyV/fkoL is rj times the work done by the driving laser field 
whose frequency is lol divided by the energy of one photon Hul [see Eq. (|17|) ]. It is a dimensionless time-dependent 
random variable, described by a probability density function P{W,T), which at least in principle can be evaluated 
based on the statistical properties of the spectral diffusion process and the stochastic Bloch equations. From Eq. 
and for a specific realization of the stochastic process Au){t), the averaged number of photons counted in time interval 
(0, T) is given by W, 

00 

(n), = ^np(n,T) = M^, (25) 

where the shot noise average is (...)s — J2^=o ■■■Pin,T). Since W is random, additional averaging over the stochastic 
process Au;(t) is necessary and statistical properties of the photon count are determined by {p(n,T)), where (• • • ) 
denotes averaging with respect to the spectral diffusion (i.e., not including the shot noise), 

/W" \ 

{p{n,T))^(—eM-W)j. (26) 



Generally the calculation of {p(n,T)) is nontrivial; however, in some cases simple behavior can be found. Assuming 
temporal fluctuations of W occur on the timescale Tc, then we have 
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(a) for counting intervals T Tc and for ergodic systems, 

lim lim i r im^iHcuL)), (27) 

which is the fluorescence hneshape of the molecule, that is, the averaged number of photon counts per unit time when 
the excitation laser frequency is lul [later we suppress lol in {I (ujl))]- Several authors fssl Fol l85>| have argued quite 



generally (though not in the SM context) that in the long measurement time limit we may use the approximation 
W ~ {I)T (i.e., neglect the fluctuations), and hence photon statistics becomes PoissonianjTOj. 

(p(n,T))^Ml):exp(-(/)T). (28) 
n] 

At least in principle (/) can be calculated based on standard lineshape theories, (e.g., in Appendix A we calculate (/) 
for our working example considered in Section 1^. Eq. (|28|) implies that a single measurement of the lineshape (i.e., 
averaged number of emitted photons as a function of laser frequency) determines the statistics of the photon count 
in the limit of long measurement time. In fact, it tells us that in this case, counting statistics beyond the average will 
not reveal any new information on the SM interacting with its dynamical environment. Mathematically this means 
that the distribution P{W, T) satisfying 

{p(n,T)) ^ dWP(W,T)—e^p(-W) (29) 
Jo ril 

converges to 

P{W,T) ^ d {W ~ {I)T) , (30) 

when T ^ oo. We argue below, however, using the central limit theorem, for cases relevant for SMS, P{W, T) is better 
described by a Gaussian distribution. The transformation Eq. H29|l is called the Poisson transform of P{W. T^^^. 
(b) in the opposite limit, T <^ Tc we may use the approximation [sg) 

W ~ loT, (31) 

where Iq — I{T — 0). In a steady state {p{n, T)) can be calculated if the distribution of intensity (i.e., photon current) 
is known 

{p{n,T))^ r dIoP{Io)e-'°^^-^^. (32) 



For example, assume that I{t) is a two state process, i. e. the case when a SM is coupled to a single slow two level 
system in a glass, then 

P{I) ^ PiS {I - h) + P26 {I - h) , (33) 

and 

(p(n, T)) - Yl P^^^ exp i-I^T) , (34) 



and if, for example, I2 ^ 0, the SM is either "on" or "off", a case encountered in several experiments |30l I32L |45 | . 
(c) a more challenging case is when T ~ Tc] later, we address this case in some detail. 

We would like to emphasize that the photon statistics we consider is classical, while the Bloch equation describing 
dynamics of the SM has quantum mechanical elements in it (i.e., the coherence). In the weak laser intensity case 
the Bloch equation approach allows a classical interpretation based on the Lorentz oscillator model as presented in 
Appendix B. 



III. SIMULATION 



To illustrate combined effects of the spectral diffusion and the shot noise on the fluorescence spectra of a SM, we 
present simulation results of spectral trails of a SM, where the fluorescence intensity of a SM is measured as a function 
of the laser frequency as the spectral diffusion proceeds 0,^3- 
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First, we present a simple algorithm for generating random fluorescence based on the theory presented in Section 
im by using the stochastic Bloch equation, Eqs. (|5|)- ((Tn|l . and the classical photon counting distribution, Eq. A 
measurement of the spectral trail is performed from t = to i = icnd- As in the experimental situation, we divide 
iend into N time bins each of which has a length of time T. For each bin time T, a random number of photon counts 
is recorded. Simulations are performed following the steps described below: 

Step (1) Generate a spectral diffusion process Aw(t) from i = to i = iond- 

Step (2) Solve the stochastic Bloch equation, Eqs. ||S||- H1U|) . for a random realization of the spectral diffusion 
process generated in Step (1) for a given value of lul during the time period < t < iond- 

Step (3) Determine W{tk,ujL) during the fcth time bin {k — 1, 2, 3, • • • N), (k — 1)T < t < kT, with a measurement 
time T according to Eq. ll^ , 

W{tk,cuL) = / dtl{t). (35) 

J(fc-1)T 

Step (4) Generate a random number < a; < 1 using a uniform random number generator, and then the random 
count n is found using the criterion, 

n— 1 n 

Y,P{j,W)<x<Y,Pij,W). (36) 
j=o j=0 



According to Eq. H23|) we find 



[n — 1)! n\ 



where r(a, z) is the incomplete gamma function. Steps (1 — 4) must be repeated many times to get good statistics. 

For an illustration purpose we choose a simple model of the spectral diffusion, which is called a two state jump 
process or a dichotomic process. We assume that the frequency modulation can be either Auj{t) — or Auj{t) = —v, 
and the flipping rate between these two frequency modulations is given by R. This model will be used as a working 
example for which an analytical solution is obtained later in Section IVI 

In Fig. n we present a simulation result of one realization of a spectral diffusion process when the fluctuation rate 
R is much smaller than F and v (slow modulation regime to be defined later). Parameters are given in the figure 
caption. Fig. Qla) shows W{tk,uJL), demonstrating the effects of the spectral diffusion process on the fluorescence 
spectra. Note that W{tk,ujL) has been defined without the shot noise. One can clearly see that the resonance frequency 
of a SM is jumping between two values as time goes on, and W shows its maximum values either at ujl — ujq = v 
or at uiL — Wo = Since shot noise is not considered in Fig. ^a), W appears smooth and regular between the 
flipping events. In Fig. db) we have taken into account the effects of shot noise as described in Step (4) and plotted 
the random counts n as a function of and t. Compared to Fig. ^a), the spectral trail shown in Fig. ^b) appears 
more fuzzy and noisy due to the shot noise effect. It looks similar to the experimentally observed spectral trails (see, 
for example, Ref. [131 )■ 

In Fig. [3 we show the evolution of the photon counting distribution P/c(n), Eq. H23|) . at a fixed laser frequency, 
chosen here as wi, — ojq — = — SF. The spectral diffusion process is identical to that shown in Fig.^ Here k denotes 
the measurement performed during the fcth time bin as described in Step (3). Notice that two distinct forms of the 
photon counting distributions appear. During the dark period at the chosen frequency — = —v in Fig. (e.g., 
2.1 X 10* <Tt < 3.5 X 10"* corresponding to 21 < < 35), Pk{n) reaches its maximum at n = with Pk{n — 0) ~ 1 
and Pfe(n > 1) 1, meaning that the probability for a SM not to emit any photon during each time bin in the dark 
period is almost one. However, during the bright period (for example, 3.5 x 10** < Ft < 4.6 x 10* corresponding to 
35 < fc < 46), Pfe(n) shows a wide distribution with {n)s — 35, meaning that on the average {n)s — 35 number of 
photons are emitted per bin during this period. As the spectral diffusion proceeds, one can see the corresponding 
changes in the photon counting distribution, Pk{n), typically among these two characteristic forms. 

In Fig. O we present a simulation result of a spectral trail for the case when the resonance frequency of a SM 
fluctuates very quickly compared with F and v (i. e. i? ^ 3> F). In this case, since the frequency modulation is 
so fast compared with the spontaneous emission rate, a large number of frequency modulations are realized during 
the time l/F, and the frequency of the SM where the maximum photon counts are observed is dynamically averaged 
between ujq — v and u>q + v (i. e. a motional narrowing phenomenon) |6^ Is^ . The width of the spectral trail is 
approximately F, and no splitting is observed even though the frequency modulation v is larger than the spontaneous 
decay rate T (v = 5T in this case). This behavior is very different from the slow modulation case shown in Fig. ^ 
where two separate trails appear at ujl = tuo ^ v- 
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FIG. 1: Spectral trails of a SM undergoing a very slow spectral diffusion process described by the two state jump model are 
shown for (a) W (without the shot noise) and (b) n (with the shot noise). Parameters are chosen as Q ^ 0.2r, R = 10"*r, 

u = 5r,T = io^r-\ tend = lo^r^^ = loor, and r = 1. 



In Fig. ^ we also show the evolution of the photon counting distribution Pk (n) during a spectral diffusion process 
at a fixed frequency ujl = loq where the lineshape reaches its maximum in the fast fluctuation case shown in Fig. |31 
Unlike the slow modulation case in Fig. |21 where one can see large fluctuations of the photon counting distributions, 
the fluctuations of the photon counting distributions are much smaller in the fast modulation case, and Pk{n) shows 
a broad Gaussian-like behavior, centered at {n)s — 12. 



IV. Q AND THREE-TIME CORRELATION FUNCTION 



Having observed the interplay between spectral diffusion and shot noise on the fluorescence spectra of a SM in 
simple simulation results of the previous section, it is natural to ask how one can analyze theoretically the photon 
counting statistics of a SM in the presence of a spectral diffusion process. The probability density of the number 
of photon counts {p{n,T)), or equivalently P{W,T) in Eq. (|29|l . would give complete information of the dynamical 
processes of a SM undergoing a spectral diffusion process, but is difficult to calculate in general. In order to obtain 
dynamical information, we will consider the mean (W) and the second moment (W'^) of the random photon counts. 

It is easy to show that the average number of photons counted in time interval (0,r) is given from Eq. H2t)|l . 



oo 

((n),)=^n(p(n,r)) = (I^), (38) 

n=0 



FIG. 2: Time evolution of photon counting distributions Pk{n) for the slow modulation case corresponding to ujl — loq — —v 
in Fig.0 Other parameters are given in Fig.0 



and the second factorial moment of the photon counts in time interval (0, T) is given by 

1)),) - l)(p(n,T)) = (W^^). (39) 

n=0 

The Mandel Q parameter is now introduced to characterize the fluctuations |70| . 



and it is straightforward to show that [7^ 



= 1 + Q, (40) 



This equation is important relating Q to the variance of the stochastic variable W. We see that Q > 0, indicating 
that photon statistics is super-Poissonian. For our classical case we anticipate: 

(a) for an ergodic system, when T/tc oo, and if Eq. Y2t^) is strictly valid, Q — > (i.e., Poissonian statistics). 
However, below we find an analytical expression for Q which is non-zero and in some cases large even in the limit of 
T —^ oo. We will discuss this subtle issue later; 

(b) in the opposite limit, T <^Tc, 

'"^ -T- (42) 



(c) if I{t) = I, independent of time, Q = 0, as expected; 

(d) it is easy to see that Q oc rj, hence when 77 — > 0, counts recorded in the measuring device tend to follow the 
Poissonian counting statistics. 
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FIG. 3: Spectral trails of a SM undergoing a fast spectral diffusion process described by the two state jump model are shown 
for (a) W (without the shot noise) and (b) n (with the shot noise). Parameters are chosen as Q — 0.2T, R — lOF, v — 5T, 
T = 10^r~\ fend = 10^r~^ = lOOT, and T = 1. 



We now consider the important limit of weak laser intensity. In this limit the Wiener-Khintchine theorem relating 
the lineshape to the one-time correlation function holds. As we shall show now, a three-time correlation function is 
the central ingredient of the theory of fluctuations of SMS in this limit. In Appendix B we perform a straightforward 
perturbation expansion with respect to the Rabi frequency ft in the Bloch equation, Eqs. (|5)l-l(Tn)l. to find 



w(i) = ^Re|^ diiexp -ij^ dt'SL{t')~^{t-ti) | 



(43) 



According to the discussion in Section^]the random number of photons absorbed in time interval (0, T) is determined 
by W = O jQv{t)dt (see Eqs. (O and (|23|)), and from Eq. lO we find 



W = — Re 
2 



t2 



j^^g-iwi,(t2-ti)-r(t2-ti)/2+i J!^ dt'Ai^it') 



(44) 



where we have neglected terms of higher order than fi^ . In standard lineshape theories Eq. H44|l is averaged over the 
stationary stochastic process and the long time limit is taken, leading to the well known result for the (unnormalized) 
lineshape 



1 -^oo 1 Z 



dre 



(45) 
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where we have set cuq = 0. The one-time correlation function C[{t) is defined by 



(46) 



where I = ±1 and (• • • ) is an average over the stochastic trajectory Aui. Eq. (|45|l is the celebrated Wiener-Khintchine 
formula relating the one-time correlation function to the average number of photon counts, i. e. the averaged lineshape 
of a SM. We now investigate lineshape fluctuation by considering the statistical properties of W. 
Using Eq. H43|l we show in Appendix B that 



exp 



diidt2dt3di4e-*'^^(*^-*^+*^-*^)-^(l*^-*^l+l*^-**l)/' 
Auj{t')dt' -i r Aw(i')di' 

Jt3 



(47) 



As can be seen from Eq. H47|l the key quantity of the theory of lineshape fluctuation is the three-time correlation 
function, 



C3 (n,T-2,r3) = (exp 



Aw(t')di' - i / Auj{t')dt' 



(48) 



which depends on the time ordering of tl,t2^ t^, 14. In Eq. H48|l we have defined the time ordered set of {ti, t2, t^, 14^} 
as {ti,tij,tiii,tiv} such that ti < tu < tm < tjv, and ti = tu — tj, T2 = tm — tu and T3 = tjv — tm. Due to 
the stationarity of the process C3 (ti, T2, r3) does not depend on the time elapsing between start of observation t = 
and ti . It has a similar mathematical structure to that of the nonlinear response function used to describe four wave 
mixing spectroscopies such as photon echo or hole burning 64J. 

In Eq. (|48|1 there are 4! = 24 options for the time ordering of (ti, t2, ^3, ^4); however, as we show below, only three 
of them (plus their complex conjugates) are needed. It is convenient to rewrite the three-time correlation function as 
a characteristic functional, 



(ri,r2,r3) = (exp 



(49) 
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m 


Sm{t) 




time ordering 


1 


1 1 

1 1 

1 

-i ' 

tl til till tiv 


P^-^(ri)P%(T2)Pl,(Ts) 


< ^2 < < t4 


2 


1 1 

' 1 

-1 1 


- T f ri 1 P^i, f r2 )P" Vrs ) 


t2 < ti < ti < 


3 


1 1 1 

U 


y -PA (n ) P'l (t2 ) Ph (t3 ) 


t2 < t\ < tz < ti 


4 


1 1 


i y p.- ) P°,, ( ra )PJT ^ (t3 ) 


t\ < t2 < ti < ts 


5 


2 1 1 

1 


lE-F»\(n)'P-i^(^2)Pw(r3) 


t2 < ts < tl < t4 


6 


-1 1 1 

-2 1 1 


hJ:P-,\r^)P^^\r2)P,-,\r,) 


ti < ti < t2 < ts 



TABLE I: Three-time correlation functions represented by six different pulse shape functions Smit) for 24 time orderings 
schemes. Values of Smit) for each time interval are shown to the left of the pulse shape. Only one representative time ordering 
scheme for each class is shown in the third column. Three other time ordering schemes belonging to the same class are obtained 
by exchanging ti <-> t4 and/or t2 t^. For example, three other time orderings belonging to m = 1 are ti < t2 < ts < ti, 
t\ < ti < t2 < ti, and < t-j, < t2 < t\. In the second column expressions of the three-time correlation functions C3" are 
shown in terms of the weight functions considered in Section|3 Note that (73""^ (n, r2, T3) and C3"(ri, r2, T3) (n=l,2,3) are 
complex conjugates of each other. 



where Sm{t) (m = 1, 2, • • • , 6) is defined in Table^as the pulse shape function corresponding to the mth time ordering. 
Let us consider as an example the case m — I, ti < t2 < < (for which ti ~ tj^i^ = t//, • • ■ , and t\ ^ t2 — t\i 
T2 = ^3 — ^2, and r3 = t4 — t^. Then the pulse shape function is given by 



-1 ti <t < tn 
Slit) = I tii<t<tui 
1 till <t< tiv, 



(50) 



and the shape of this pulse is shown in the first line of Table Similarly, other pulse shapes describe the other time 
orderings. 

The four dimensional integration in Eq. (|47|l is over 24 time orderings. We note, however, that 



dt'Aw(t') - i / dt'Atj(t') 



g -lWi(i2-tl+t3-t4)-r(|tl-t2| + |t.3-t4|)/2 ^ 



in Eq. H47() has two important properties: (a) the expression is invariant when we replace ti with ^4 and t2 with ^3, 
and (b) the replacement of ti with t2 (or ^3) and of ^3 with ti{ov ti) has a meaning of taking the complex conjugate. 
Hence it is easy to see that only three types of time orderings (plus their complex conjugates) must be considered. 
Each time ordering corresponds to different pulse shape function, Sm{t)- In Tabled for all six time ordering schemes, 
the corresponding pulse shape functions are presented. We also give expressions of C™(i) for the working example to 
be considered soon in Section Hvl 
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1 h(t) 



FIG. 5: A schematic representation of the spectral diffusion process modeled by two state jump process(random telegraph 
noise). The fluctuating transition frequency is given by Aaj(t) = vh{t). 



We note that if ri = T3 the pulse in Eq. (|5()|l is identical to that in the three-time photon echo experiments. The 
important relation between lineshape fluctuations and nonlinear spectroscopy has been pointed out by Plakhotnik in 
Ref. |68| in the context of intensity-time-frequency-correlation measurement technique. 



V. TWO STATE JUMP MODEL: EXACT SOLUTION 



A. Model 



In order to investigate basic properties of lineshape fluctuations we consider a simple situation. We assume that 
there is only one bath molecule that is coupled to a chromophore by setting J — 1 and Auji{t) — vh{t) in Eq. Q), 
where v is the magnitude of the frequency shift and measures the interaction of the chromophore with the bath, and 
h{t) describes a random telegraph process h(t) = 1 or h{t) = —1 depending on the bath state, up or down, respectively. 
For simplicity, the transition rate from up to down and vice versa are assumed to be R. The generalization to the 
case of different up and down transition rates is important, but not considered here. A schematic representation of 
the spectral diffusion process is given in Fig. 

This model introduced by Kubo and Anderson in the context of stochastic lineshape theory [60I l6ll | is called the 
sudden jump model, and it describes a stochastic process that describes fluctuation phenomena arising from Markovian 
transitions between discrete states jo^, 0|- For several decades the Kubo- Anderson sudden jump model has been a 
useful tool for understanding lineshape phenomena, namely, of the average number of counts (n) per measurement time 
T, and has found many applications mostly in ensemble measurements, e.g., NMRI6.2'] and nonlinear spectroscopy '645 . 
More recently, it was applied to model SMS in low temperature glass systems in order to describe the static properties 
of lineshapes fM ITU ITa | and also to model the time-dependent fluctuations of SMS|6^|63,E1|- In this paper, we will 
consider this model as a working example to study properties of lineshape fluctuations. 

The above model can describe a single molecule coupled to a single two level system in low temperature glasses 
as explained in Section FVIII In this case v depends on the distance between the SM and the two level svstem|88|. 
Another physical example of this model is the following: consider a chromophore that is attached to a macromolecule, 
and assume that conformational fluctuations exist between two conformations of the macromolecule. Depending on 
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the conformation of the macromolecule, the transition frequency of the chromophore is either luq — v oi cuo + i^|31 



B. Solution 



By using a method of Suarez and Silbey|89(, developed in the context of photon echo experiments, we now analyze 
the properties of the three-time correlation function. We first define the weight functions, 



exp 



"'0 



(51) 



where the initial (final) state of the stochastic process Aui is i (/) and a = or a = ±1 or a = ±2. For example, 
P++(t) is the value of {e^^o '^'^(* )'^* ) for a path restricted to have Atj(O) = v and Aivit) — v. The one-time correlation 
function G\{t) defined in Eq. (|46|l can be written as sum of these weights. 



(52) 



where a prefactor i is due to the symmetric initial condition and the summation is over all the possible paths during 
time t{\. e. i = ±, j = ±). Also by using the Markovian property of the process, we can express all the functions 
in terms of the weights. For example, for the pulse shape in Eq. H5U|) 



(ri, r2, ra) = J E (^i) P% (^2) Pli (rg) , 



(53) 



ij.k.l 



where the summations are over all possible values of i — ±, j — ±, fc = ±, and I — ±. The other C™ functions are 
expressed in terms of weights in Table ^ Explicit expressions of the weights for the working model are given by 



P|_(t) = Fl+(t) = i?exp(-i?t)sin(Yii)/Yi 

sin {Yity 



Pl±{t) = exp(-i?t) 



cos (Yit) =F iiy- 



Yi 



(54) 
(55) 
(56) 



Pi'/it) = C.C. Plf (<) where C.C. denotes complex conjugate, and Yi = Vi^^ - R^- P^'^{t) is given by the same 



expressions as the corresponding Pfi^{t) with v replaced by 2v. 



Now we can evaluate {W) and {W"^) explicitly. First we consider {W"^) and, in particular, focus on the case m = 1, 
< h < t2 < h < ti < T. By using Table H the contribution of {W'^) to Eq. (HTJ is 



— / du I dt 











3 / dt2 dtie~ 
"'0 



Jt(t2-tl+t3-*4)-r(|tl-t2| + lt3-U|)/2 



i,j,k,l 

We use the convolution theorem of Laplace transform four times and find 



1234 



+ zc.i)Pj],(.s)P,i (s + r/2-z^i) i. 



where C ^ denotes the inverse Laplace transform, where the Laplace T —> s transform is defined by 



/(.'^) = / fiT)e 
Jo 



-"'^dT, 



(57) 



(58) 



(59) 
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and the Laplace transforms of the functions P^j{s) are listed in Eqs. (|133|I - H137() . Considering the other 23 time 
orderings we find 

(ir^) = ^ 3 V I fc' (, + r/2 + i^j p' («) p*' (, + r/2 - i^^) 
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+Pij' is + r/2 + z^l) Pji. (s + r) p+i (,s + r/2 - z^^) 

+P±i (s + r/2 - ic^i) Pji (s) P+i (s + r/2 - tujL) 

(s + r/2 - iu;l) Pji (s + r) P+1 (.s + r/2 - i^l) 
+2P+1 (s + r/2 - ilul) P+k' [s + T- 2iujL) P^i^ [s + V/2 - ilol) + C.C.] } 



(60) 



Eq. H6Q(I . which can be used to describe the lincshapc fluctuations, is our main result so far. In Appendix A we invert 
this equation from the Laplace s domain to the time T domain using straightforward complex analysis. Our goal is 
to investigate Mandel's Q parameter, Eq. (|41|l : it is calculated using Eqs. (|60|l and 



{W) = ^£-1 \^^Y^P7.^s + T/2 + iujl) + C.C. j , (61) 

which is also evaluated in Appendix A. As mentioned, Eq. H61|l is the celebrated Wiener-Khintchine formula for the 
lineshape (in the limit of T — > oo) while Eq. H6U|) describes the fluctuations of the lineshape within linear response 
theory. Note that {W'^) and {W) in Eqs. and (|FT|l are time-dependent, and these time-dependences are of interest 
only when the dynamics of the environment is slow (see more details below). Exact time-dependent results of (W"^) 
and (W), and thus Q, are obtained in Appendix A. The limit T oo has, of course, special interest since it is used in 
standard lineshape theories, and does not depend on an assumption of whether the frequency modulations are slow. 
The exact expression for Q in the limit of T ^ cx3 is given in Appendix B, Eqs. (|150|l and (|151() . These equations are 
one of the main results of this paper. It turns out that Q is not a simple function of the model parameters; however, 
as we show below, in certain limits, simple behaviors are found. 



VI. ANALYSIS OF EXACT SOLUTION 



In this section, we investigate the behavior of Q{u>l,T) for several physically important cases. In the two state 
model considered in Section in addition to two control parameters ujl and T we have three model parameters that 
depend on the chromophore and the bath : E, R, and i^. Depending on their relative magnitudes we consider six 
different limiting cases: 

1. i?< r < 1/ 

2. i?< i/< r 

3. r<i/<i? 

4. i/<r<i? 

5. r<i?<iy 

6. i^<i?<r 

We discuss all these limits in this section. 



A. Slow Modulation Regime : R <t. iy,r 



We first consider the slow modulation regime, R <ti where the bath fluctuation process is very slow compared 
with the radiative decay rate and the frequency fluctuation amplitude. In this case the foregoing Eqs. (|60ll and 161(1 
can be simplifled. This case is similar to situations in several SM experiments in condensed phases 6, 7, 8, 9, 12, 13]. 

Within the slow modulation regime, we can have two distinct behaviors of the lineshape depending on the magnitude 
of the frequency modulation, v. When the frequency modulation is slow but strong such that R <^ T <^ i/ [case 1], 
the lineshape exhibits a splitting with the two peaks centered at = zLu (see Fig. (S^a)). On the other hand, when 
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the frequency modulation is slow and weak, i? ^ <C F [case 2] a single peak centered at = appears in the 
lineshape (see Fig. [TJa)). From now on we will term the case F ^ as the strong modulation limit and the other 
case F ^ J/ as the weak modulation limit. The same distinction can also be applied to the fast modulation regime 
considered later. 

In the slow modulation regime, we can find Q(T) using random walk theory and compare it with the exact result 
obtained in Appendix A. In this regime, the molecule can be found either in the up + or in the down — state, if 
the transition times (i.e., typically 0{l/R)) between these two states are long; the rate of photon emission in these 
two states is determined by the stationary solution of time-independent Bloch equation in the limit of weak laser 
intensitylz^ [see also Eqs. and (jHH)] 



I±{^l)^ _ , r., . (62) 



Now the stochastic variable W = I{t)dt must be considered, where I{t) follows a two state process, I{t) = /+ or 

I{t) = I_ with transitions H > — and > + described by the rate R. One can map this problem onto a simple two 

state random walk problem 99| , where a particle moves with a "velocity" either or /_ and the "coordinate" of the 
particle is W. Then from the random walk theory it is easy to see that for long times (T 3> ^/R), 

{W) ^ [^^^) T = {I)T, (63) 

meaning that the line is composed of two Lorcntzians centered at ±z/ with a width determined by the lifetime of the 
molecule, 

{I) = lil++I-), (64) 
and the "mean square displacement" is given by 



{W') - ^ J'^ T. (65) 
After straightforward algebra Eqs. (|63|) . H65|) . and (|41|) yield 

in the limit T ^ 1/R. The full time-dependent behaviors of (W) and (VF^) are calculated in Appendix C using the 
two state random walk model,,9Qj. From Eq. (|161|l . we have Q as a function of the measurement time in the slow 
modulation regime 



2 



where the "variance" of the lineshape is defined by 

((A/)^> = ^[(/+ + (68) 

Eq. (|67|l shows that Q is factorized as a product of frequency dependent part and the time-dependent part in the 
slow modulation regime. It is important to note that Eq. 1)67(1 could also have been derived from the exact result 
presented in Appendix A by considering the slow modulation conditions. Briefly, from the exact expressions of (W) 
and (VF^) in the Laplace domain given in Appendix A, by keeping only the pole s's that satisfies |Re(s)| ~ 0{TZ) and 
neglecting other poles such that |Re(s)| S> R, we can recover the result given in Eq. l|H7|l. thus confirming the validity 
of the random walk model in the slow modulation regime. In the limits of short and long times, we have 

ii^R-^ T » 1/R 



(I) 



(69) 

jj^T T « 1/R . 

We recover Eq. JHEl) from Eq. in the limit of T ^ oo. Eq. for T < 1/i? is a special case of Eq. Note 
that the results in this subsection can be easily generalized to the case of strong external fields using Eqs. H96|) and 
(El. 
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FIG. 6: Case 1 {R <ti T <^ i/) in the steady-state limit. Exact results for lineshape, Eq. I|95|l . and for Q, Eqs. H15()|l . are 
shown for the strong, slow modulation case as functions of ujl- Parameters are chosen as v = IGHz, F — lOOMHz, Q = F/10, 
R = IHz — lOOHz, and T — > oo. They mimic a SM coupled to a single slow two level system in a glass. Note that lineshape 
does not change with R while Q does. 



1. strong modulation : case 1 

When the T ^ ly, case 1, the two intensities I± are well separated by the amplitude of the frequency modulation, 
2i>, therefore we can approximate (/) ~ /+/2 when /_ ~ and vice versa, which leads to 

((A/)2) ~ (/)2. r«i. (70) 

In this case Q is given by 

f {I)R-^ T> 1/i? 

(71) 

[ {I)T T < 1/i? . 

We show the lineshapes and Q as functions of the laser frequency at the steady-state limit (T oo) for case I 
in Fig. |S| In all the calculations shown in figures of the present work, we have assumed an ideal measurement, r] — 1. 
Values of parameters are given in the figure caption. This is relevant for the case that a chromophore is strongly 
coupled to a single two level system in a low temperature glass ||S,'^i|. Since Q ~ in this case, both the lineshape 

(Fig. EJa)) and Q (Fig. E^b)) are similar to each other; two Lorentzian peaks located at uj^ = ±;y with widths F. 
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Note that in this hmit the value of Q is very large compared with that in the fast modulation regime considered later. 
While the lineshape is independent of R in this limit, the magnitude of Q decreases with i?, hence it is Q not (/) 
which yields information on the dynamics of the environment. 



2. weak modulation : case 2 



On the other hand, in case 2, where the fluctuation is very weak {v <^V), we notice that /+ ~ /_ from Eq. H62|l 
(note that when there is no spectral diffusion, = 0, /+ = thus (5 = 0). The lineshape is given by a single 
Lorentzian centered at = 0, 

(/) - (72) 



Since ^ F in this case we expand I± in terms of v to find 

'd{I)y _ ly^n 
d^J ~4(^+r2/4)4 



.9. 9 I d{I)\ ly^fl'^T^ivf , , 

((A/)2) ^ I _W j ^«r (73) 



T^(^)Vi T»l/R 

, (74) 

T r<i/i? . 



Therefore, in case 2, Q is given by 



Note that Q oc ■^(^^)^ in the weak, slow modulation case while Q oc (/) in the strong, slow modulation case. 

In Fig. 13 we show the lineshape and Q for the weak, slow modulation limit, case 2. The lineshape (Fig. |7| (a)) 
is a Lorentzian with a width F, to a good approximation, thus the features of the lineshape do not depend on the 
properties of the coupling of the SM to a environment such as ly and R. On the other hand, Q (Fig. [7| (b)), shows a 
richer behavior. Recalling Eq. (|74|l for T ^ 1/R in Eq. (|73|1 . it exhibits doublet peaks separated by ~ C(F) with the 
dip located aX lol = 0, and its magnitude is proportional to 1/R. We will later show that this kind of a doublet and 
a dip in Q is a generic feature of the weak coupling case, found not only in the slow but also in the fast modulation 
case considered in the next section. Note that when the SM is not coupled to the environment (i^ = 0), Q = as 
expected. 

Both in the strong (Fig. 0) and the weak (Fig. [T)) cases, we find Q{lol = 0) ~ 0. This is expected in the slow 
modulation case considered here. Physically, when the laser detuning frequency is exactly in the middle of two 
frequency shifts, -izv, the rate of photon emissions is identical whether the molecule is in the up state or in the 
down state {—v). Therefore, the effect of bath fluctuation on the photon counting statistics is negligible, which leads 
to Poissonian counting statistics at cj^ = 0. 



3. time-dependence 



Additional information on the environmental fluctuations can be gained by measuring the time dependence of Q 
in the slow modulation regime. In Fig. |S|we show Q versus the measurement time T for the strong, slow modulation 
limit, case 1, both for the exact result calculated in Appendix A and for the approximate result, Eq. (|67|) . We choose 
the resonance condition ojl — v-, and used parameters relevant to SMS in glass systems. The approximate result 
based on the two state random walk model, Eq. (|57|) . shows a perfect agreement with the exact result in Fig.|HIa) as 
expected. When T <C 1|R^ Q increases linearly with T as predicted from Eq. (|69|l for RT <^ 1, and it reaches the 
steady-state value given by Eq. Hti9|l when T becomes RT ^ 1. We also notice that even in the long measurement 
time limit the value of Q is large: Q = 'b x lO^ry in the example given in Fig.|S^a) (including the detection efficiency). 
Therefore, even if we consider the imperfect detection of the photon counting device (for example, 77 — 15% has been 
reported recently 91]), deviation of the photon statistics from Poissonian is observable in the slow modulation regime 
of SMS. This is seemingly contrary to propositions made in the literature ^3,, 33., 85j, in which the claim is made 
that the Poissonian distribution is achieved in the long time limit. We defer discussion of this issue to the end of this 
section. We note that it has been shown that Q{T) can be very large (Q ~ iC) at long times in the atomic three-level 
system with a metastable state but without a spectral diffusion process [s^- Fig.|S{b) illustrates that the steady-state 
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FIG. 7: Case 2 {R <ti i/ ^ T) in the steady-state limit. Exact results for lineshape, Eq. H95|l . and for Q, Eqs. H15()|l . are 
shown for the weak, slow modulation case as functions of ujl- Parameters are chosen as i/ = IMHz, F = lOOMHz, fl — F/IO, 
i? = IHz - lOOHz, and T ex. 



value of Q{T) is reached when RT ^ 1, and the magnitude of Q in the steady-state decreases as 1/i? as predicted 
from Eq. HH9|I for RT ^ 1. This therefore illustrates that valuable information on the bath fluctuation timescales can 
be obtained by measuring the fluctuation of the lineshape as a function of measurement time. 

Fig. El shows a two-dimensional plot of Q as a function of the laser frequency and the measurement time for the 
parameters chosen in Fig. |HIa). We observe that two Lorentzian peaks located at resonance frequencies become 
noticeable when T ~ 1/i?. Similar time dependent behavior can be also found for the weak, slow modulation limit, 
case 2, however, then Q along the uil axis shows doublet peaks separated by ~ 0(r) as shown in Fig. [3 

We can extend our result and describe photon statistics beyond the second moment. Based on the central limit 
theorem, the probability density function of the two state random walk variable W in the long time limit is described 
by 



P{W,T) ^-^^ G{W,T) 



1 



^/4:TTDT 



exp 



{w ~ vry 



(75) 



with D ((A/)2)/(2i?) and V (/) when 77 = 1. We also note that Q ~ 2D/V in the long time hmit. By using 
Eq. ^ 



{pin,T)) 



dWG{W,T)—eyip{-W). 



(76) 
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FIG. 8: Case 1 (7? ^ F <g z/) in the time-dependent regime, (a) Q{T) is shown both for the exact result and for the 
approximate result given in the linear-log scale. This figure demonstrates that our exact results given in Appendix A are well 
reproduced by the approximation, Eq. H67|l . The same parameters are used as in Fig. Elexcept ul =v = IGIfz, R = IHz, and 
T (varied), (b) Q(r) for different values of R in the slow modulation regime in the log- log scale. All the parameters except for 
R are chosen the same as those in (a). 
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FIG. 9: Case 1 {R <Si T <^ v) in the time dependent case. Two dimensional plot of Q{ujl,T) for the slow modulation as a 
function of and T. The same parameters are used as in Fig. |5Ja) except lul (varied). 

Eq. H76() shows that in the long time limit the photon statistics is the Poisson transform [6^ of a Gaussian. This 
transformation can be found explicitly (one can slightly improve this approximation by considering a normalized 
truncated Gaussian with G{W, T) — for W < 0; we expect that our approximation will work well when V^T ^ D). 
In contrast to Eq. H75|l . the proposed Eq. H3U|I suggests that G{W,T) be replaced with a delta function, for which a 
single parameter V controls the photon statistics, while according to our approach both V and D, or equivalently, (W) 
and Q are important. Mathematically, when T —^ oo the Gaussian distribution, Eq. H75I) may be said to "converge" 
to the delta function distribution, Eq. (|30|) . in a sense that the mean is (W) ~ 0{T) while the standard deviation is 
AW = {{W^) — (W)"^)^^^ ^ OiT^/"^). This argument corresponds to the proposition made in the literature, Eq. l|5n)l. 
and then the photon statistics is only determined by the mean {W) ~ VT . However, physically it is more informative 
and meaningful to consider not only the mean lyV) but also the variance (AVF)^ ~ 2DT ~ Q{W) since it contains 
information on the bath fluctuation. In a sense the delta function approximation might be misleading since it implies 
that Q = in the long time limit, which is clearly not true in general. 

B. Fast Modulation Regime : !^,r < i? 

In this section we consider the fast modulation regime, v,T <^ R. Usually in this fast modulation regime, the 
dynamics of the bath (here modeled with R) is so fast that only the long time limit of our solution should be 
considered [i.e. Eqs. H15U|) and H151|) ]. Hence all our results below are derived in the limit of T — > oo, since the time 
dependence of Q is irrelevant. The fast modulation regime considered in this section includes case 3 (F <C <C i?) as 
the strong, fast modulation case and case 4 (z^ ^ F ^ i?) as the weak, fast modulation case. 

It is well known that the lineshape is Lorentzian^4j in the fast modulation regime (soon to be defined precisely). 



(77) 
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where Tf is the Hne width due to the dephasing induced by the bath fluctuations and given by 

T/ = ^. (78) 

The hneshape given in Eq. IjTTII exhibits motional narrowing, namely as R is increased the line becomes narrower and 
in the limit R ~^ oo the width of the line is simply F. 

Now we define the fast modulation regime considered in this work. When i? — > cxd (and other parameters fixed) 
Q ^ 0, so fluctuations become Poissonian. This is expected and in a sense trivial because the molecule cannot 
respond to the very fast bath, so the two limits i? — > oo and v = (i.e., no interaction with the bath) are equivalent. 
It is physically more interesting to consider the case that R oo with Tf/T remaining finite, which is the standard 
definition of the fast modulation regimej^. In this fast modulation regime the well known lineshape is Lorentzian 
as given in Eq. l(77|) with a width F + F/. 

Here we present a simple calculation of Q in the fast modulation regime based on physically motivated approxi- 
mations. We justify our approximations by comparing the resulting expression with our exact result, Eqs. I|15UI) and 
(|151|l . When the dynamics of the bath is very fast, the correlation between the state of the molecule during one pulse 
interval with that of the following pulse interval is not significant in the stochastic averaging Eq. (|49|) . Therefore 
in the fast modulation regime, we can approximately factorize the three-time correlation function C™(Ti,r2,r3) in 
Eq. I|49|) in terms of one-time correlation functions C{ (r) as was done by Mukamel and Loring in the context of four 
wave mixing spectroscopy l92| , 

^~ia jf^" Ai^{t')dt'-ib f^^'^" Auj{t')dt' -ic f*^'^^ Aw(t')dt'^ 
g-ia ///^ A^(t')dt'\^ (^^-^bf,\'^" A^(t')dt'\^ ^e~^' 

= Cnri)Cl{T2)Ct{Ts), (79) 

where a, 6, and c are values of Smit) in time intervals, tj < t < tjj, tu < t < tju, and tjjj < t < tjy, respectively. 
For example, a = —1, 6 = 0, and c — 1 for the case m = 1 (see Tabled). The one-time correlation function is 
evaluated using the second order cumulant approximation, which is also valid for the fast modulation regime |9^. 

C[{t)~expl-l^ f dt' f dt" {Alu it' ) Auj {t" ))] . (80) 



\ Jo Jo ) 

Note that these two approximations, Eqs. (|79|l and H80() . are not limited to the two state jump model but are generally 
valid in the fast modulation regime. The frequency correlation function is given by 

(At^(t')Aw(t")> - exp(-2i?|t' - t"|) (81) 
for the two state jump model. Substituting Eq. (|81() into Eq. 1801) we have 

t e-2«*-r" 



C\(t) ~ exp 



(82) 



2R {2Rf 

In the limit of i? — )■ oo we can approximate C\(f) further 

C{(i) ~ exp(-;2ryt/2), (83) 

and in turn evaluate all the three-time correlation functions corresponding to different time orderings within these 
approximations using Eqs. (|79|l and (|83|l . We note that the lineshape shows a Lorentzian behavior with the width 
F + F/ in Eq. H77|l because of the exponential decay of the one-time correlation function in Eq. H83|) . Since we have 
simple analytical expressions of the three-time correlation functions, it is straightforward to calculate Q in the fast 
modulation regime. After straightforward algebra described in Appendix D, Q in the fast modulation regime is 

^^Vf \%Yf^l + 2(r + F^)(8r2 + 15FFy + SF^^i + F/(F + 2F/)(F + V ff 

Q = ^— 



(r+2r/ 



(84) 
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Note that the exact expression for Q obtained in Appendix B yields the same expression as Eq. (|84(l when i? — > oo 
but T f remains finite, which justifies the approximations introduced. 

Let us now estimate the magnitude of these fluctuation. We consider lol = Q since the photon current is strongest 
for this case (i.e., the lineshape has a maximum at cj^ = 0), then 

Q = , (85) 

r(r + r^f (r + 2r^) 

where the detection efficiency 77 has been restored. The maximum of Q is found when F/ = r(l + •\/5)/2 and then 
Q ~ Q.iQXrjVl^ /T^ . Even if we take \^\/T = 1/10 and = 5 x 10^^ as reasonable estimates for a weak laser field and 
detection efficiency we find Q ~ 2 x 10^"' ^ 1 which shows the difhculties of measuring deviations from Poissonian 
statistics in this limit. We note, however, that values of Q as small as \Q\ < 10"'' have been measured recently 
(although in the short time regime, not in the steady-state case)|55. Therefore it might also be possible to observe 
the deviation from Poissonian photon statistics in the fast modulation limit under appropriate experimental situations. 

In Fig. ^] we show the results of exact steady-state calculations for the lineshape and Q [Eq. (|150|l and (|151|l ] 
for different values of the fluctuation rate R in the fast modulation regime. We have chosen the parameters as 
1/ = lOOMHz, F = IMHz, and fl = F/10, and R is varied from IGHz to lOOGHz, corresponding to case 3. For this 
parameter set we have checked that the fast modulation approximation for Q given in Eq. (|84|l agrees well with the 
exact calculation, Eqs. (|150() . (|151() . 

The lineshape shown in Fig. llOf a) shows the well known motional narrowing behavior. When R = IGHz, the 
lineshape is a broad Lorentzian with the width F-l-F/ — F/ (note that Fy ^ F in this case). As R is increased further 
the line becomes narrowed and finally its width is given by F as i? — > 00. 

Compared with the lineshape, Q in Fig. llOf b) shows richer behavior. The most obvious feature is that when R —> 00, 
Q 0. This is expected since when the bath is very fast the molecule cannot respond to it, hence fluctuations are 
Poissonian and Q ^ 0. It is noticeable that, unlike Q shows a type of narrowing behavior with splitting as R 

is increased. The lineshape remains Lorentzian regardless of Fy in the fast modulation case, while Q changes from a 
broad Lorentzian line with the width T f when F/ ^ F to doublet peaks separated approximately by F when F/ <C F, 
that will be analyzed in the following [see, for example, Eq. (|86|l ]. Therefore, although the value of Q is small in 
the fast modulation regime, it yields additional information on the relative contributions of F and Fy which are not 
differentiated in the lineshape measurement. 



1. strong modulation : case 3 

We further analyze the case that the bath fluctuation is both strong and fast, case 3, F <C ^ i?. The results 
shown in Fig. ^| correspond to this case. We find that Eq. H84I) is further simplified in two different limits, F/ ^ F 
and F/ < F, 

When F/ ^ F, both Q and the lineshape are a Lorentzian located at cj^ = with a width F/, which yields the relation 
Q ~ 2{I{ujl))/^, and both exhibit motional narrowing behaviors. In the other limit, F/ ^ F, we have neglected 
©((Fy/F)^) terms with additional conditions for utj^, ^ <C ^ T~- ^^^^^ '^^^'^ *3 shows a splitting behavior 

at \ujl\ ^ C(F). We note that Eq. 1)86(1 for F/ ^ F is the same as Eq. (174(1 This is the case because the very fast 
frequency modulation corresponds to the weak modulation case, if we recall that the dephasing rate due to the bath 
fluctuation is given by F/ = jR in the fast modulation regime. 

In Fig. we have checked the validity of the limiting expressions of Q for the Lorentzian and the splitting cases 
(Eq. ((86(1 by comparing to the exact results, Eqs. ((150() . ((151() . In Fig. Ilir a). the parameters are chosen such that 
F/ = lOOF while in Fig. llir bl F/ = F/100. Approximate expressions (dashed line) show a good agreement with exact 
expressions (solid line) in each case. 



2. weak modulation : case 4 



Now we consider the weak, fast modulation case, case 4 (z/ <C F ^ i?). In this limit, the lineshape is simply a single 
Lorentzian peak given by Eq. (1771) with F -I- F/ ~ F. We obtain the following limiting expression for Q from Eq. (184(1 . 
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FIG. 10: Case 3 {F <ti v <^ R) in the steady-state limit. Lineshape and Q in the last modulation regime are shown as lunctions 
ol LJL. Parameters are chosen a.s u = lOOMHz, F = IMHz, Q = T/IO, R = IGHz - lOOGHz, and T ~* oo. 



noting Tf = i^'^/R, 



^. (87) 



Also here Q shows splitting behavior. It is given by the same expression as that in the strong, fast modulation case 
with F/ <C r (see Eq. (|86|l ) and also as that in the weak, slow modulation case with T ^ 1/R (see Eq. (|74|l ) (see also 
Table CTl in Section IVIPp . When Q is plotted as a function of uil, it would look similar to Fig. irTT b). 



C. Intermediate Modulation Regime 

So far we have considered four limiting cases: (i) strong, slow, (ii) weak, slow, (iii) strong, fast, and (iv) weak, fast 
case. We now consider case 5 {T <^ R <^ v) and case 6 (:^ <C i? <C F). They are neither in the slow nor in the fast 
modulation regime according to our definition. 
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FIG. 11: Case 3 {T <^ u <^ R) in the steady-state limit. Exact calculations (solid line) of Q, Eqs. 1150^ and 1)151^ are compared 
with the approximate expressions, Eq. 18611 . (dashed line) for (a) the Lorentzian (F/ ^ F), and (b) the sphtting (F/ <^ F) cases 
in the fast modulation case. Parameters are chosen as v = IGHz, F = IMHz, Q = F/10, T — > oo, and R = lOGHz in (a) while 
R= lO^GHz in (b). 



1. case 5 



In case 5(r <C i? ^ J^), the bath fluctuation is fast compared with the radiative decay rate but not compared with 
the fluctuation ampHtude. Because T R, i/ in this case we can approximate the exact results for (/) and Q by their 
limiting expressions corresponding to F ^ 0, yielding Eqs. H100() and l|152(l . and an important relation holds in this 
limit, 

lim Q = - lim(/). (88) 

Note that the same relation between Q and (/) was also found to be valid in one of the fast modulation regimes, 
Eq. 186|l with F/ 3> F. Fig. 1121 shows that in this case the limiting expressions approximate well the exact results. 



2. case 6 



In case 6, since ^ i? <C F we can approximate the exact results by considering small v limit in Eqs. I|15()(l and 
(|151|l for Q, and Eq. H95|) for the lineshape. By taking this limit, we flnd that the lineshape is well described by a 
single Lorentzian given by Eq. (|72|) . and Q by Eq. (|74|l . We note that for all weak modulation cases, cases 2, 4, and 6, 
the lineshape and Q behave in a unique way described by Eq. l|7^ and Eq. l(7^ for RT » 1, respectively, and both 
the slow and fast modulation approximate results are valid in this case (see Table lll)l . Also a simple relation between 
(J) and Q holds in the limit, v ^ 

limQ^^limi-f^V. (89) 



^0 R v^Q (!) \duj 



The exact results of lineshape and Q in Fig. E| show good agreement with the approximate results, Eqs. ((7^ and 
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FIG. 12: Case 5 [T <^ R u) in the steady-state limit. Exact results for lineshape, Eq. I|95|l . and for Q, Eqs. H15()|l . p51^ . 
are compared with approximations corresponding to F — > hmit, Eqs. I|100|l and l|152|l . respectively. Parameters are chosen as 
= IGHz, F = 5MHz, Q = F/10, R = lOOMHz, and T ^ oo. 



D. Phase Diagram 

We investigate the overall effect of the bath fluctuation on the photon statistics for the steady-state case as the 
fluctuation rate R is varied from slow to fast modulation regime. To characterize the overall fluctuation behavior of 
the photon statistics, we define an order parameter q, 

g = — / Q{ul)Awl, (90) 

where Q in the steady state is given in Eqs. H150() and H151() . Before we discuss the behavior of q it is worthwhile 
mentioning that the lineshape is normalized to a constant regardless of R and F, 

/oo ttO^ 
^(/K))dc.i = — , (91) 

which can be easily verified from Eq. H45|l . In contrast, q exhibits nontrivial behavior reminiscent of a phase transition. 
In Figs. 1141 and 1151 we show q versus R/T. The figures clearly demonstrate how the photon statistics of SMS in the 
presence of the spectral diffusion becomes Poissonian as — > oo or ^ (i.e. q ^ when i? ^ oo or ^ 0). 
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FIG. 13: Case 6 {v <ti R <^ T) in the steady-state limit. Exact results for lineshape, Eq. I|95|l . and for Q, Eqs. H15()|l . p51^ . 
are compared with approximations corresponding to — > hmit, Eqs. H72^ and l|89|l . respectively, Parameters are chosen as 
= IMHz, r = lOOMHz, n = r/10, R = IOMHz, and T ^ oo. 



We first discuss the strong modulation regime, v ^T, shown in Fig. 1 141 with v/T = 100 in this case. Depending on 
the fluctuation rate, there are three distinct regimes: 

(a) In the slow modulation regime, R -^T, q decreases as 1/R. The approximate calculation (dot-dashed line) based 
on the slow modulation approximation, Eq. H71|l shows good agreement with the exact calculation. 

(b) When R is such that F <^ R ly [case 5], the intermediate regime is achieved, and q starts to show a plateau 
behavior. The plateau behavior is found whenever Q = 2{I{uj]^)) /T, which yields g — 1 as can be easily seen from 
Eqs. H9UI) and H91|l. As R is increased further such that v ^ R, the fast modulation regime is reached, and q still 
shows a plateau until R ~ iP /T. In this regime the Lorentzian behavior of Q is observed in Eq. H8()|l for F/ F. 

(c) When the bath fluctuation becomes extremely fast such that R ^ i^^/F, the splitting behavior of Q is observed, 
as discussed in Eq. I|8f>(l for F/ <C F, and then q oc 1/R similar to the slow modulation regime. The approximate value 
of q based on fast modulation approximation, Eq. H84|l (dotted line) shows good agreement with the exact calculation 
found using Eqs. (|150|l and (|151|) . Finally, when R oo, Q = 0. As mentioned this is expected since the molecule 
cannot interact with a very fast bath hence the photon statistics becomes Poissonian. 

Now we discuss the effect of magnitude of frequency fluctuation, i/, on q in Fig. ^1 We have calculated g as a 
function of 1/R as is varied from I'/T — 100 (circles) to z//F = 0.01 (downward triangles). We see in Fig. 1151 that 
the (g, i?/F) diagram exhibits a behavior similar to a phase transition as ly is varied. In the strong modulation regime 
(j//F ^ 1), three distinct regimes appear in the (g,i?/F) diagram {R <T, T < R< f^/T, and R > /V), while in 
the weak modulation regime {v/T <C 1) g always decreases as 1/R. When three parameters, F, and i?, have similar 
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FIG. 14: q versus R/T in the strong modulation case. Filled circles are the result of exact calculation based on Eqs. and 
when iz/T = 100. The dot-dashed and dashed curves are calculations based on approximate expressions of Q, Eq. H71|l 
(slow modulation case) and Eq. (I84II (fast modulation case), respectively. 



values, J/ ~ r ^ i?, there appears a "critical point" in the "phase diagram". 

Table ^ summarizes various expressions for the lineshape and Q found in the limiting cases of the two state jump 
model investigated in this work. For simplicity, we have set 77 = 1. We see that although the fluctuation model 
itself is a simple one, rich behaviors are found. We believe that these behaviors are generic (although we do not 
have a mathematical proof). In the weak modulation cases, both (/) and Q can be described by a single expression, 
irrespective of the fluctuation rate R. However, in the strong modulation cases, (/) and Q change their qualitative 
features as R changes. 



VII. CONNECTION TO EXPERIMENTS 

Single molecule spectroscopy has begun to reveal the microscopic nature of low temperature glasses [ol IT3 . ITsl IT^ 
l65l | . An important question is whether the standard tunneling model of low temperature glass developed by Anderson, 
Halperin, Varma|93| and Phillips "9^ is valid or not. As far as macroscopic measurements of acoustic, thermal, and 
optical properties are concerned, this model has proved to be compatible with experimental results 88]. However, on 
a more microscopic level we do not have much experimental or theoretical proof (or disproof) that the model is valid. 
At the heart of the standard tunneling model is the concept of a two level system(TLS). At very low temperatures, 
the complicated multidimensional potential energy surface of the glass system is presumed to reduce to a multitude 
of non-interacting double well potentials whose two local minima correspond to reorientations of clusters of atoms or 
molecules [93. Hence the complicated behavior of glasses is reduced to a simple picture of many non-identical a,nd 
non-interacting TLSs. For a different perspective on the nature of low temperature excitation in glasses, see Ref. ,9^. 

Geva and Skinner [l^ have provided a theoretical interpretation of the static lineshape properties in a glass (i.e. (W)). 
The theory relied on the standard tunneling model and the Kubo-Anderson approach as means to quantify the 
lineshape behavior (i.e. the time dependent fluctuations of W are neg lected). In Ref. O, the distribution of static 
lineshapes in a glass was found analytically and the relation of this problem to Levy statistics was demonstrated. 

Orrit and coworkers |8j ^] have measured spectral trails as well as the lineshapes and the fluorescent intensity 
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FIG. 15: Phase diagram of q versus R/T. Symbols are the results of exact calculations based on Eqs. 115UII and 115111 as u/T is 
varied from the lower-left to the upper-right cases: 0.01 (downward triangle), 0.1 (upward triangle), 1 (diamond), 10 (square), 
100 (circle). The dashed curves are calculations based on appropriate approximate expressions of Q for each case, Eqs. IHZJ 
(weak modulation), 1711 (strong, slow modulation), and 1841 1 (strong, fast modulation). 



correlation function g^'^\t). In Ref. |Q spectral trails of 70 molecules were investigated and 22 exhibited behaviors 
that seemed incompatible with the standard tunneling model. While the number of molecules investigated is not 
sufficient to determine whether the standard tunneling model is valid, the experiments are approaching a direct 
verification of this model. 

Our theory in the slow modulation limit describes SMS experiments in glasses. The TLSs in the glass flip between 
their up and down states with a rate R determined by the coupling of TLS to phonons, the energy asymmetry, and 
the tunneling matrix element of the TLS|U0|. When a TLS makes a transition from the up to the down state, or 
vice versa, a frequency shift u oc occurs in the absorption frequency of a SM, where r is the distance between 
the SM and the TLS. The dependence is due to an elastic dipole interaction between the SM and the TLS. In a 
low temperature glass, the density of TLSs is very low, hence one finds in experiment that the SM is coupled strongly 
to only a few TLSs. In some cases, when one TLS is in the vicinity of the SM, it is a reasonable approximation to 
neglect all the background TLSs. In this case, our theory describes SMS for chromophores in glasses with a single 
TLS strongly coupled to SM. Extension of our work to coupling of SM to many TLSs is important, and can be done 
in a straightforward manner provided that the TLSs are not interacting with each other. 

Fleury et al.^ measured g'''^\t) for a single terrylene molecule coupled to single TLS in polyethylene matrix. They 
showed that their experimental results are well described by 



{R+1+ 

where Rj^ and R- are the upward and downward transition rates, respectively. These two rates, R+ and i?_ arc due 
to the asymmetry of the TLS. This result is compatible with our result for Q in the slow modulation limit. When 
Eq. (|92|l is used in Eq. (0 with i?+ = i?_ = R for the symmetric transition case considered in this work, we exactly 
reproduce the result of Q{T) for the slow modulation given in Eq. (|67|l . [The asymmetric rate case i?+ ^ i?_ is also 
readily formulated for Q in the slow modulation limit, and again leads to a result compatible with Eq. H92|l .] Hence, 
at least in this limit, our results are in an agreement with the experiment. 
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intermediate 


fast 


weak 


i? < < r : case 2 


«^ < -R < r : case 6 


!^ < r < i? : case 4 




1 




2 


strong 


-R < r < i/ : case 1 


r < i? < 1/ : case 5 


r < < i? : case 3 


/n ^+ + ^- 




, n'(r + r/) 















TABLE II: Parameter regimes investigated in this work are classified with expressions for the lineshape, (/), and Q for each 
case. I± have been defined in Eq. 1)62^ . 



As far as we are aware, however, measurements of photon counting statistics for SMS in fast modulation regimes 
have not been made yet. The theory presented here suggests that even in the fast modulation regime, the deviation 
from Poisson statistics due to a spectral diffusion process might be observed under suitable experimental situations, 
for example, when the contributions of other mechanisms such as the quantum mechanical anti-bunching process and 
the blinking process due to the triplet state are known a priori. In this case, it can give more information on the 
distribution of the fluctuation rates, the strength of the chromophorc-cnvironmcnt interaction, and the bath dynamics 
than the lineshape measurement alone. 

VIII. FURTHER DISCUSSIONS 

All along in this work we have specified the conditions under which the present model is valid. Here we will 
emphasize the validity and the physical limitations of our model. We will also discuss other possible approaches to 
the problem at hand. 

In the present work, we have used classical photon counting statistics in the weak laser field limit. In the case of 
strong laser intensity, quantum mechanical effects on the photon counting statistics are expected to be important. 
From theories developed to describe two level atoms interacting with a photon field in the absence of environmen- 
tal fluctuations, it is known that, for strong field cases, deviations from classical Poissonian statistics can become 
significant [toI l97| . One of the well known quantum mechanical effects on the counting statistics is the photon anti- 
bunching effect |55l 

lIlEaizllzlIzi In this sub-Poissonian behavior is obtained, — 1 < Qqm < 0, where 

the subscript "qm" stands for quantum mechanical contribution. It is clear that when the spectral diffusion process 
is significant (i. e. cases 1 and 2) any quantum mechanical correction to Q is negligibly small. However, in the fast 
modulation case where we typically found a small value of Q, for example, Q ~ 10~^, quantum mechanical correc- 
tions due to anti-bunching phenomenon might be important unless experiments under extremely weak fields can be 
performed such that |Qqm| < Qsd, where the subscript "sd" stands for spectral diffusion. The interplay between truly 
quantum mechanical effects and the fast dynamics of the bath is left for future work. 

In this context it is worthwhile to recall the quantum jump approach developed in the quantum optics community. 
In this approach, an emission of a photon corresponds to a quantum jump from the excited state to the ground state. 
For a molecule with two levels, this means that right after each emission event, pee = (i. e. the system is in the 
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ground state). Within the classical approach this type of wavefunction collapse never occurs. Instead, the emission 
event is described with the probability of emission per unit time being Tpeeit), where Peeit) is described by stochastic 
Bloch equation. At least in princip le, the quantum jump approach, also known as the Monte Carlo wavefunction 
approach "o?, ^QSl lo^ llOOl llOll Il02j| . can be adapted to calculate the photon statistics of a SM in the presence of 
spectral diffusion. 

Another important source of fluctuation in SMS is due to the triplet state dynamics. Indeed, one of our basic 
assumptions was the description of the electronic transition of the molecule in terms of a two state model. Blinking 
behavior is found in many SMS experiments J.8, 19, 20, 25, 26]. Due to the existence of metastable triplet states 
(usually long lived) the molecule switches from the bright to the dark states (i. e. when molecule is shelved in the 
triplet state, no fluorescence is recorded). Kim and Knight jsO] pointed out that Q can become very large in the case 
of the metastable three level system in the absen ce of spectral diffusion. This is especially the case when the lifetime 
of the metastable state is long. Molski et aZ. |53 . IT03 | have considered the effect of the triplet state blinking on the 
photon counting statistics of SMS. 

Therefore at least three sources of fluctuations can contribute to the measured value of Q in SMS ; (i) Qam, well 
investigated in quantum optics community, (ii) Qtripict, which can be described using the approach of Ref. 80], and 
now we have calculated the third contribution to Q, (iii) Qsd- Our approach is designed to describe a situation for 
which the spectral diffusion process is dominant over the others. 

It is interesting to see if one can experimentally distinguish Qsd and Qtripict for a SM in a condensed environment. 
One may think of the following gedanken experiment: consider a case where the SM jumps between the bright to the 
dark state, and assume that we can identify the dark state when the SM is in the metastable triplet state. Further, 
let us assume that dark state is long lived compared to the time between emission events in the bright state. Then, 
at least in principle, one may filter out the effect of the dark triplet state on Q especially when the timescale of the 
spectral diffusion process is short compared with the dark period by measuring the photon statistics during the bright 
period. 



IX. CONCLUDING REMARKS 



In this paper we have developed a stochastic theory of single molecule fluorescence spectroscopy. Fluctuations 
described by Q are evaluated in terms of a three-time correlation function C3 {ti,T2,T3) related to the response 
function in nonlinear spectroscopy. This function depends on the characteristics of the spectral diffusion process. 
Important time ordering properties of the three-time correlation function were investigated here in detail. Since the 
fluctuations (i.e., Q) depend on the three-time correlation function, necessarily they contain more information than 
the lineshape which depends on the one-time correlation function Ci (ti) via the Wiener-Khintchine theorem. 

We have evaluated the three-time correlation function and Q for a stochastic model for a bath with an arbitrary 
timescale. The exact results for Q permit a better understanding of the non-Poissonian photon statistics of the single 
molecule induced by spectral diffusion. Depending on the bath timescale, different time orderings contribute to the 
lineshape fluctuations and in the fast modulation regime all time orderings contribute. The theory predicts that Q is 
small in the fast modulation regime, increasing as the timescale of the bath (i.e. 1/-R) is increased. We have found 
nontrivial behavior of Q as the bath fluctuation becomes fast. Results obtained in this work are applicable to the 
experiment in the slow to the intermediate modulation regime (provided that detection efhciency is high), and our 
results in the more challenging fast modulation regime give the theoretical limitations on the measurement accuracy 
needed to detect Q > 0. 

The model system considered in this work is simple enough to allow an exact solution, but still complicated enough 
to exhibit nontrivial behavior. Extensions of the present work is certainly possible in several important aspects. It 
is worthwhile to consider photon counting statistics for a more complicated chromophore-bath model, for example, 
the case of many TLSs coupled to the chromophore, to see to what extent the results obtained in this work would 
remain generic. Also the effects of a triplet bottleneck state on the photon counting statistics can be investigated as a 
generalization of the theory presented here. Another direction for the extension of the present theory is to formulate 
the theory of SMS starting from the microscopic model of the bath dynamics (e. g. the harmonic oscillator bath 
model) . Effects of the interplay between the bath fluctuation and the quantum mechanical photon statistics on SMS 
is also left for future work. 

The standard assumption of Markovian processes (e.g. the Poissonian Kubo-Anderson processes considered here) 
fails to explain the statistical properties of emission for certain single "molecular" systems such as quantum dotsl 2lll2a. 
|23 ]. Instead of the usual Poissonian processes, a power-law process has been found in those systems. For such highly 
non-Markovian dynamics stationarity is never reached and hence o ur a pproach as well as the Wiener-Khintchine 
theorem does not apply. This problem has been investigated in Ref. |108| . 
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X. APPENDIX A: CALCULATION OF LINESHAPE 



In this appendix, we calculate the lineshape for our working example. We set J = 1 and luq ^ and as mentioned 
the stochastic frequency modulation follows Auj{t) = vh{t), where h{t) describes a two state telegraph process with 
h(t) = +1 (up) or h{t) = — 1 (down). Transitio ns from s tate up to down and down to up are described by the rate R. 

We use the marginal averaging aDProach |8ll Il04 Il05l | to calculate the average lineshape. Briefly, the method gives 
a general prescription for the calculation of averages {y) {y is a vector) described by stochastic equation y = M( t)y, 
where M(t) is a matrix whose elements are fluctuating according to a Poissonian process (see Refs. |^ Iin5l | for 
details). We define the marginal averages, {v{t))x , {u{t))x and {'w{t))x where x = + or x ^ — denotes the state of 
the two state process at time t. For the stochastic Bloch equation, the evolution equation for the marginal averages 
are 
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(93) 



where St = SlT v. The steady state solution is found by using a symbolic program such as Mathem,atica \li)(\ . 

{I {lol)) = \n {{v,t)- + {vst) + ) - {{wst)+ + (wst)- + 1) , (94) 
where {{wst)+ + (wst)- + 1) /2 represents the steady state occupation of the excited level. We find 

{I{cuL)) = i, (95) 



A = rn^ [(r + 2R) {t{4ujI + (r + 4Rf) + Av^{t + aR)) + 2n'^v{v + aR)] , 

B = T{T + 2R) ((4(wi - v^) - r(r + AR)f + 16^1 (F + 2R)) + A^f {AujIT{T + 3i?) 

+(r2 + AYR + 4t/2)(r2 + 3ri? + ar^)) + aq.^t{t + ar). 

Remark 1 When i? — > 0, it is easy to show that (/ (wl)) is a sum of two Lorentzians centered at ±1^, 

{I{lol)) = ^{v++v^), 
where v± are steady state solution of Bloch equation for two level atom (see Ref. [z^]) 

17 r/2 



v± 



2 (,5±)2+r2/4 + r!V2' 



(96) 



(97) 



Remark 2 If — > 0, the solution can also be found based on the Wiener-Khintchine formula, using the weights 
Pj^^(s) defined in Appendix A. 



(/(c.O) = ^Re 



^ Pr^'{iujL+Tl2) 



(98) 
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which gives 



{A{wI-u^)-T{T + AR)Y + IQloI[T + 2RY- ^ ' 

Now if r ^ we get the weh known result of Kubo 

= (100) 

then in the slow modulation limit, v ^ R the line (/ {lol)) exhibits splitting (i.e., two peaks at ±i/) while for the fast 
modulation regime v <^ R the line is a Lorentzian centered at oj^ =0 and motional narrowing is observed. 
Remark 3 If ^ the solution reduces to the well known Bloch equation solution of a stable two level atom, which 
is independent of R. 

Remark 4 We have assumed that occupation of state + and state — are equal. More general case, limited to weak 
laser intensity regime, was considered in Ref. pflf. 

XI. APPENDIX B: PERTURBATION EXPANSION 

In this appendix, we find expressions for the photon current using the perturbation expansion. We also use the 
Lorentz oscillator model to derive similar results based on a classical picture. 

We use the stochastic Bloch eauation |8ll Is^ to investigate Pee{t) in the limit of weak external laser field when we 
expect Pee — 0, Pgg ~ 1 for times TT > 1. We rewrite Eqs. (P)- ((TI7|l 

= -Tpee + Y {peg - Pge) , (101) 
^ - -[Z5i(t)+r/2]pg,-y (pee-Pss), (102) 

where peg = Pe^e*"^*, pge = PgeC"*"^', pee = Pee and pgg = Pgg. Using Eqs. UlUlfl . (11021) . the normahzation condition 
pee + pgg = 1, and pge = C.C[peg\, the four matrix elements of the density matrix can be determined in principle 
when the initial conditions and the stochastic trajectory Aw(t) are specified. We use the perturbation expansion 

Peeit) = pi?(0 + ^P^^^t) + f^2p(2)(t) + . . . , (103) 
PAt) = PfKt) + ^P^gKt) + ^'Pfe\t) + ■■■: (104) 

and initially p^e}{t = 0) = pi"^{t = 0) = for i = 1, 2, 3, • • • . We insert Eqs. ((TU5|l . (Pljl into Eqs. pUT|) .pU21). and 
first consider only the zeroth order terms in fl. We find pW {t) = /?ee(0) exp(— Ft), this is expected since when the 
laser field is absent, population in the excited state is decreasing due to spontaneous emission. The off-diagonal term 
is p^g} it) — /Oge''(0) exp[— ri/2 — i |C ! At'SjXt')], this term is described by the dynamics of a Kubo- Anderson classical 
oscillator, x = [-i5L{t) - r/2]a;j63t. 

First Order Terms 

The first order term is described by the equation 



dt 



This equations yields the solution 



One can show that p^ee [t) is unimportant for times Vt ^ 1, like all the other terms which depend on the initial 
condition. 



Plin0 = e-" / e"'dt'Imb(°)(Ol. (106) 
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For the ofF-diagonal term we find 



dt 



z<54i)+r/2]p«-^(2p(°)-l) 



The transient term plPj — pee(O) exp(— Fi) is unimportant, and Eq. (|107|l yields 



P« it) - \ 



At I exp 



(107) 



(108) 



Using v{t) = mm p^g} {t) we find Eq. 

According to the discussion in the text the number of photons absorbed in time interval (0, T) is determined by 

time integration of the photon current ~ dtv{t). Using Eq. (|43|l and definition of 5L{t), we obtain Eq. 144|l . 
It is convenient to rewrite Eq. (|44|) also in the following form 



4 Jo 



^^^g-iwi,(t2-ti)-r|t2-ii|/2+i//^2 Acj(t')dt' 



(109) 



We calculate the average number of counts (W) using Eq. 144|l . The integration variables are changed to t — t2~ti 
and ti, and for such a transformation the Jacobian is unity. Integration of ti is carried out from to T — r, resulting 
in 



(W) = —Re 



|^dr(l-^)e— ^-/^Cr^(r) 



(110) 



and C[{t) = l^e-^iio ^^(t'W) is the one-time correlation function. In the limit of T ^ oo we find Eq. (|45|) . 
Using Eq. H109(l the fluctuations are determined by 



16 Jo 



T pT pT rT 



X < exp 



^0 JQ 

i 



diidi2dt3di4e-^"^(*--'i+*^-*3)-r(l*i-*2|+l*3-t4|)/2 
dt'Auj{t')+i [ dt'Auj{t') 



(111) 



and changing integration variables, ^3 and ^ t^, yields Eq. (|47|l . We replaced ^3 and to get a pulse shape 

similar to that in the three-time photon echo experiment (when ii < t2 < ^3 < ti)- Notice that the derivation did 
not assume a specific type of random process ALu{t) and our results are not limited to the two state telegraph process 
we analyze in the text. 



Second Order Terms 



According to Eq. |(T7J| and ^ F Pee{t)dt = Q v{t)dt > 1. We now show that this equation is valid within 
perturbation theory. For this aim we must consider second order perturbation theory. The equation for the second 
order term 



dpi'e'(t) 



dt 



-r/5, 



n(2) 



Im 



P'^jit) 



yields 



-r*Re 



\\t2 ^i^^^iii+t2)/2-iiJL{t2-ti)+iS'^ Auj(t')dt' 

Jo 

-)2~(2), 



(112) 



(113) 



where we have neglected terms depending on initial condition. Since Pedt) — fl pee (t) we see that the response 
is quadratic with respect to the Rabi frequency as we expect from symmetry (population in excited state does not 
depend on sign of Eq). Using Eq. (|113|l we find 



Pee{t) = 2^ e 



2-rt 



dt2 
Jo 



(J^^gr(ti+t2)/2-iWi(t2-tl)+i/t'= Au{t'}dt' 



(114) 
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('2') 

Now, we consider the standard ensemble measurement and average pie yt) with respect to history of the process 
Aa;(t'). Assuming stationarity and a changing variables, we find in the limit of T — > oo 



{piV (oo)) = —Re 



dTCr^(r)e-^^/2-^"^^ 



(115) 



which is the lineshape. From Eq. H45|l we see that r(/9ie'(oo)) = limT-.oo(VF/r) = (/(wl)). Thus the theory of the 
averaged lineshape can be based on either v (first order perturbation theory) or on p^e (second order perturbation 
theory) . 

Instead of the averages, let us consider the stochastic variables Ffi^ pl?J {t)dt and fl^ v^^\t)dt, where w*^^^ = 



Im 



p^eg [t) ■ Using Eq. (|113|l . and the Laplace T ~* s transform 
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/•oo 

/sto is) = / die-^*Re 
Jo 

is a functional of the stochastic function Auj{t). Using Eq. (|44|l we find 
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(116) 



(117) 



(118) 



Comparing Eq. (|116|l and Eq. (|118|l we see that a theory based on p^e or on v are not entirely identical. However, for 
long times FT :s> 1, we may use small s behavior (i.e., F/ (s + F) ~ 1) and F J^^ peeit)dt — Q v{t)dt as expected. 



Classical Lorentz Model 



The stochastic Bloch equation is a semiphenomenological equation with some elements of quantum mechanics in 
it. To understand better whether our results are quantum mechanical in origin we analyze a classical model. Lorentz 
invented the theory of classical, linear interaction of light with matter. Here we investigate a stochastic Lorentz 
oscillator model. We follow Allen and Eberley[107j who considered the deterministic model in detail. The classical 
model is also helpful because its physical interpretation is clear. We show that for weak laser intensity, the stochastic 
Bloch equations are equivalent to classical Lorentz approach. 

We consider the equation for harmonic dipole |e|a;(i) in the driving field, £{t) — E cos{ujLt), 

lei 

i + Tx + Lu^{t)x ^ —£{t) (119) 
m 

where all symbols have their usual meanings and ujo{t) — loq + vh{t) is a stochastic time-dependent frequency and 
|i^ft.(t)| ^ LJo- AH along this section we use symbols which appear also in the Bloch formalism since their meanings in 
the Bloch and in the Lorentz models are identical, as we show below. 
We decompose x into two parts. 

x{t) — xq [u cos {ujLt) — V sin {uJLt)] , (120) 

xq is a time-independent constant, while u and v vary slowly in time. The work done by the laser force F = 
\e\E cos{u)Lt) on the particle is dW — ~Fdx hence ^ — —Fx. Using Eq. (|120|l we find 

dW 1 

— - -\e\ExoUJLV (121) 

and the overbar denotes average over rapid laser oscillations (e.g., we assume that the noise term h{t) evolves slowly 
if compared with the laser period 2'k/ujl). Since h{t) is stochastic so is v{t) hence the power is also a stochastic 
function. 



36 



As in Ref. |107| . we assume that u and v vary slowly in time such that 



\u\<t:ujL\vl |u|<wi|u|, \v\ <t: ujl\u\, \v\<^Lui\v\, 
then insert Eq. H12U|) into Eq. (|119|) and find two equations for the envelopes u and v, 

u = [ujL- ujo[t)\ V ~ —, 



[UJL - LOoit)] U - 



Tv 



\e\E 



(122) 

(123) 
(124) 



where the relation [wo(^) — '^i]/(2^L) — '^o(^) — was used. Comparing Eqs. (|121|l . (|123() . H124I) with Eqs. l(TB|l . 
|H|l- (ll()ll . we will now show that in the weak laser intensity limit the Bloch equation describes the dynamics described 
by the Lorentz model. To see this clearly, note that when — > 0, pee — 0, and hence w ~ —1/2. Therefore if we 
replace —ilw in the Bloch equation, Eqs. (|S|)- (|10|I with 0/2, the Bloch equations for u and v become uncoupled from 
that for w. Using this approximation we find 



= -5L{t)u 



Tu 
T' 



(125) 
(126) 



It is clearly seen that the Bloch equation in the weak intensity limit [Eq. H125(l and H126|l ] has the same structure 
as the Lorentz equation [Eq. H123|) and H124|) ]. Note that two parameters, xq and m only appear in the Lorentz model 
while two other parameters, 17 = — deg • E/fi, and K in the Bloch equation. To make the equivalence between these two 
approaches complete, the following relations can be deduced by comparing Eqs. (|123l) . H124|) with Eqs. (|125ll . (|126|l . 
and Eq. ifHT]! with Eq. lfTH)l. 



mujLXo ' 



■f-^ — -|e|a;o, fko^ -ma;|a;g. 



To conclude, when the laser intensity is not strong the stochastic phenomenological Bloch equation describes classical 
behavior. 



XII. APPENDIX A: EXACT CALCULATIONS OF (W) AND (W^) 
In this appendix we use straightforward complex analysis and find 



where 



a (5 
Us 

Us 



i,j,k.l 
i,j.k.l 

^ E P^ + ^i^ + r) Pki + S-) , 



i,j,k,l 



= - E + -'-^Pjk i-s + 2.S-) Pli (s + s_) 



(127) 

(128) 
(129) 
(130) 
(131) 
(132) 



ijj^k^l 
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and s± = r/2 ± iuiL. The Laplace transforms of the weights P^jir) are 



P^_{s 



P\{s 

PiL{s 

p-lis 



R + s 
s (s + 2R) ' 
R 

s (s + 2R) ' 
R + s — iiy 

(S-Sl)(s-S2)' 

pi + (,s)=P_,-^(,s) 
i? + s + w 

(s - Si) (s - S2)' 



Plj {s)\iy^2iy, 



(133) 
(134) 
(135) 
(136) 

(137) 



where 



Sl,2 



-R± yfW^. 



The inverse Laplace transforms of ^i(s) in Eq. (|132ll are calculated using standard methods of complex analysis to 
yield e,(T): 



m— 1 

where the simple poles zi are given by 
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(138) 



R-y/W^, 



and /i(2:, r) is 



T 1 - e 



zT 



/i(z,T) = ----.- 



(139) 



2z z2 

Notice that if ujl = 0, zi = and Z2 = 2:4, then the poles become second order. Also, we can neglect exponential 
decays exp(zir) for i = 1 — A since FT ^ 1, and the term exp(z5T) is important when RT < 1. 



UT)=Y1 [s'l2(s)(s-z™) 



h{Zni,T), 



where zi — — s_|_ + si, Z2 = + ,S2, 2:3 = — s_ + si, 24 = — s_ + -2:5 = — F, = — 2i? — F and 

T 1 - e^^ 
/2(2,T) = -i-^^. 

The expression 5^^2(5) (s — Zm) l^^^ in Eq. (|140(l is the residue of s'^^2{s) when s = z^. 



(140) 



(141) 



UT) = [s3(s-Z3)6(s)l /l(-3,T) 
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E [^'(^-^™)'C3(S) 
m— 1 
2 



E 



ds 



s^(s-z^)2^3(s) 



.hizTn,T) 



fl{Zm,T), 



(142) 
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where zi = — ,s_ + si, Z2 = — s_ + S2, zj, = —2R and 



2T Te^^ 3(1 - e^^) 
/3(.,T) = ^ + -^ + -^ + 



(143) 



2 



S — Zr, 
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rn—l 

4 



d 

dl 



S^(S - Z™)^^4(s) 



m— 3 



f2iZm,T) 

i — Zj-n 

f2{Zm,T), 



where zi = — s_ + si, Z2 = — s- + S2, 23 = — F, Z4 = —2R — T and 



, , , T Te^^ 2(1 - e^^) 



Finally, 



?5(T) = ^ [s2(s-z™f|5(a 



m — 1 
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m—1 

4 



ds 



S^(S - Zm)^C5(s) 



m— 3 



/4(^,r) 

i — Zj-ri 



(144) 



(145) 



(146) 



where Zi = — s_ + si, Z2 = + ^2; -23 = — 2s_ + si and Z4 — — 2s_ + S2 a-nd si.2 = — -R ± y/B? — Ai'^ . 
From Eq. (|110|l the average counting number can be written as 



(W) 



£-1 



1 



C^\s + s+) + C.C. 



4 i^s 

which leads to Eq. (|61|l with Eq. (|52|l . It is also easy to show that 



{w)^—c- 

4 I Sl - S2 



1 



[-S2/2 {si - s+, T) + S1/2 (S2 - s+, T)] + C.C. 



where f2{z, T) was defined in Eq. H141() . In the limit of large T we have 

^/^7(s+) + c.c. 



rj2 



(147) 



(148) 



(149) 



XIII. APPENDIX B: Q IN THE LONG TIME LIMIT 

The exact expression for the Q parameter in the long time limit is given by Q = Numerator[(5]/Denominator[(3], 
where 

Denominator [Q] ^ TR [F^ + BF^i? + IQv^R + 4F {v^ + + ul)] 

X [F^ + 4F^i? + 16Fi? + ujI) + 4F2 {iv"^ + R^ + 2uol) 
+16 {v^ - 2iy^Lul + R^Lol + luI)] 

X [F'' + 8F^i? + 32Fi? ^ ^2 ^ 8p2 f^^2 _^ + uol) 

+l&{v^ ^2v^ul+AR^uol+ul)]^ (150) 
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and 



Numerator[Q] = dii^^n^ [r"t^i + 28r^°Rujl 

+16r®i? {U6R^ujI + 21luI + ly^ {AR^ + 21ujI)) 
+hl2T'^v^R^ {llv^ + Q2R^ujI + mR^ul + 3w^ 
+1/"* (I8i?2 - 19w|) + 5z/2 {im^ul + cj^)) 
+ 16r^ (620i?^wi + 183i?^a;^ + Qui 
W (4i?2 + &ujI) + t/2 (25i?'' + l82R^Lul + 4wi)) 
+32r'^i?('832i?''tj? +424i?2t^f + 42tj? 



'1/ 
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(25i?2 + 42w|) + (38i?^ + 437i?2w|, + 28^^)) 

'4 D ,6 / r)'? d2 I 1,1,, 2 \ 1 o/in d2 , .2\ 



-128T^R (t/^ (27i?2 + Uivl) + 2 {lOR^ + 7ujI 



X {AR'^LUL + 2iy^ (28i?^ + 65R'^luI - IujI) 



(8i?^ + 625R'^ujI + 207R^Loi - 14a;|)) 
-1024rz/2i?2 (•j.s _^ 8i?4t^i + 3i?2wi + ujI 
fi/^ (7i?2 - AujI) + [-WR^ujI + 6cji) 
f (40i?*t^i + i?2wi - 4tji)) 
-64r5 (688i?^w|, + 5527?"^! + IWR^Ji + 4^1 

(28i?^ + 653i?*wi + \Q2R^u:\ - 4w^)) 

-2048j/2i?3 [yi _ 4j,6^2 _^ 4^4^4 ^ 5^2^6 _^ ^8 

Fi/"* {hR^uj\ + Gcj^) - 2t/2 {^R^uj\ + 2w^)) 
-256r3 (^i/*^ (4i?2 + ) + t^2 |-^2 _^ ^2 |-4^2 ^ ^ 
fi/*^ (39i?* + 4i?2t^2 _ ^ 
Fi/^ (20i?^ + WAR^uPi - IR^ujI + 6cji) 



-V 



(356i?^w| + 215i?''wi - UR^ujI - 4w|))] . (151) 



To obtain these results wc find the small s expansion of the Laplace transforms (W"^) and {W). These give in a 
standard way the long time behavior of the averages (W''^) and {W) with which Q is found. These equations were 
derived using Mathematica (lQ(3\ . without which the calculation is cumbersome. Note that Q is an even function of ul 
and V, as expected from symmetry. For some special cases, discussed in the text, the exact results are much simplified. 
Remark 1 In the limit of F — > 0, which corresponds to the limit of Kubo's lineshape theory, we get 



2i/2r!2i? 



r[(tj2 _j,2)2_^4^2^2] 



(152) 



Remark 2 For cj^ = 0, we get 



g ^ 2,m^.^RiT + 2R) 

F (F + 4i?) (r2 + 4i/2 + 2Fi?) (F2 + 4j/2 + 4ri?)2 



XIV. APPENDIX C: Q IN THE SLOW MODULATION REGIME 



In the two-state random walk model where the "velocity" of the particle is either or /_ , we can conveniently 
calculate the first and second moments of the "coordinate" W by introducing the generating function G(k,T), 




(154) 
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where k is an auxiliary variable. For the stochastic process where the "velocity" of the particle alternates between 
between and /_ with transition rate R, we can easily evaluate G{k, T) using the Laplace s — > T transformation, 



Now from Eq. (| 15411 wc have (W) and {W^) 

{W) = 

By taking derivatives of G(k,T) we find 

dG{k,T) 



s + ik{I+ + I^)/2 + 2R 
{s + ikl+ + i?)(s + ikl_ +R)- 



dG{k,T) 



dik 



k=0 



/(<)dt 



d^G{k,T) 



d{iky 



k=0 



dik 
d^Gik,T) 



£-1 



fe=0 



which yield after the inverse Laplace transform, 



(W) = 



t2 . t2 



s2(s + 2i?) s^{s + 2R) 



(/++/- )2 ^2 , fT e-2«T 



1 



R 



(155) 

(156) 
(157) 

(158) 
(159) 

(160) 
(161) 



Q is given by Eq. (|67|l in the slow modulation regime. 

XV. APPENDIX D: Q IN THE FAST MODULATION REGIME 

Based on the approximation introduced in the text we can calculate Q in the fast modulation regime. Once 
the factorization of the three-time correlation functions is made in Eq. (|79|l . ^i(s), the functions determining (W"^) 
[Eq. (|132|l ] can be written as 
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+ s+)C°{s)Cl{s + s-), 
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+ s+)CUs + T)Cl{s + s^), 
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s^)Cl{s + 2s^)Cl{s + s^), 


(166) 



where s± = T/2 ± iuL as defined in the Appendix A. Then {W) and {W'^) are calculated from Eqs. and ifnTjl . 

{W) = ^/:-i|^i^i(.)|, (167) 
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^-^ {>(.)} 



(168) 
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where 

Fi{s) = ci(s) + C.C., (169) 



F2{s) = (^l + --^j^i(s)^ + 2s(ci(s)^C2(s) + C.C.), (170) 

ci(s) = Ci(s + s_), (171) 
C2(s) EE Cf{s + 2s_). (172) 

After the cumulant approxhnation is made for C{(t) and the long time hmit is taken in Eq. H83(l ci(z) and €2(2:) are 
simply given by 

- . + ._Vr,/2 ' (1^3) 
^ . + 2.\2r/ (1^4) 

Since only the long time limit is relevant for the calculation of Q in the fast modulation regime we make expansions 
of (W) and {W'^) around s = 0, and find in the long time limit 

n' Fm-Fiio)F{io) 
^ = T fM ■ ^ ^ 

Note that Eq. (|175|l is valid once the factorization approximation is made irrespective of the second order cumulant 
approximation. After performing a lengthy but straightforward algebra using Eqs. H172|l - H175(l . we obtain the result 
of Q in the fast modulation regime given in Eq. (fH^ . 
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